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Chapter 1 

Gravity and Geometry 



1.1 The gravitational force 

Gravity is the most universal force in nature. As far as we can tell from observations and 
experiments every object, every particle in the universe attracts any other one by a force 
proportional to its mass. For slow moving bodies at large distances this is a central force, 
inversely proportional to the square of the distance. As the action is reciprocal, and since 
according to Newton action and reaction forces are equal in magnitude, the expression for 
the gravitational force between two objects of mass Mi and M2 at a distance R is then 
determined to have the unique form 



The constant of proportionality, Newton's constant of gravity, has dimensions of accelera- 
tion per unit of mass times an area. Therefore its numerical value obviously depends on 
the choice of units. In the MKS system this is 



It is also possible, and sometimes convenient, to fix the unit of mass in such a way that 
Newtons constant has the numerical value G = 1. In the natural system of units, in which 
also the velocity of light and Planck's constant are unity: c = h = 1, this unit of mass is 
the Planck mass mp: 



apart and do not move too fast with respect to one another. In particular, it describes 
the motions of celestial bodies like the moon circling the earth, or the planets orbiting the 
sun, as well as those of terrestrial objects like apples falling from a tree, or canon balls 



G = 6.672 59(85) x lO^^^m' 



kg ^ s ^. 



(1.2) 




1 



2 



in free flight. Ever since Newton tliis unification of celestial and terrestial mechanics has 
continued to impress people and has had a tremendous impact on our view of the universe. 
It is the origin and basis for the belief in the general validity of physical laws independent 
of time and place. 

1.2 Fields 

Although as a force gravity is universal, Newton's law ( |1.1| ) itself has only limited validity. 
Like Coulomb's law for the electrostatic force between two fixed charges, Newton's law 
holds strictly speaking only for static forces between bodies at rest. Moreover, and unlike 
the electric forces, there are modifications at smaller, finite distances which can be observed 
experimentally. 

For example, if the gravitational force would have a pure 1/ dependence, the orbits of 
particles around a very heavy central body would be conic sections: ellipses, parabola's or 
hyperbola's, depending on the energy and angular momentum, in accordance with Kepler's 
laws. The observation of an excess in the precession of the perihelion of the orbit of Mercury 
around the sun by LeVerrier in 1845, and improved by Newcomb in 1882 was one of 
the first clear indications that this is actually not the case, and that the gravitational force 
is more complicated. 

The exact form of the gravitational forces exerted by moving bodies is a problem with 
many similarities to the analogous problem in electrodynamics. The understanding of 
electro dynamical phenomena greatly improved with the introduction of the concept of 
local field of force. This concept refers to the following characteristics of electrodynamical 
forces: 

• the influence of electric charges and currents, and of magnetic poles, extends through- 
out empty space; 

• the force on a standard test charge or test magnet at any given time depends on its 
location with respect to the source of the field and the relative state of motion; 

• changes in the sources of the fields give rise to changes in the force on test objects at 
a later time, depending on the distance; the speed of propagation of disturbances in 
empty space is finite. 

In the case of electrodynamics this speed turned out to be a universal constant, the speed 
of light c. One of the most striking consequences of these properties, which follow directly 
from the mathematical description of the fields as expressed by Maxwell's equations, is the 
existence of electromagnetic waves, local variations in the fields which transport energy, 
momentum and angular momentum over large distances through empty space. Maxwell's 
predictions were magnificently verified in the experiments of Hertz, confirming the reality 
and characteristics of electromagnetic waves |@]. From these results it became clear, for 
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example, that light is just a special variety of such waves. Therefore optical and electro- 
magnetic phenomena have a common origin and can be described in a single theoretical 
framework. 

The concept of field has found other applications in physics, for example in the phe- 
nomenological description of fluids and gases. Treating them as continua, one can describe 
the local pressure and temperature as fields, with local variations in these quantities prop- 
agating through the medium as sound waves or heat waves. 

As in electrodynamics, also in gravity the concept of field has taken a central place. 
What Maxwell achieved for electromagnetism, was accomplished by Einstein in the case 
of gravity [Q: to obtain a complete description of the forces in terms of space- and time- 
dependent fields. This field theory is known as the general theory of relativity. One of the 
most remarkable aspects of this theory is, that it provides an interpretation of gravitational 
phenomena in terms of the geometry of space-time. Many detailed presentations of the 
geometrical description of space-time can be found in the literature, including a number of 
books included in the references at the end of these lecture notes 0-|]14|. A general outline 
of the structure of the theory is presented in the following sections. 



1.3 Geometrical interpretation of gravity 

The first step in describing a local field, both in electrodynamics and gravity, is to specify 
the potentials. The force on a test particle at a certain location is then computed from 
the gradients of these potentials. In the case of gravity these quantities have not only a 
dynamical interpretation, but also a geometrical one. Whereas the electromagnetic field is 
described in terms of one scalar and one vector potential, together forming the components 
of a four vectorQ y4^(x), the gravitational field is described by 10 potentials, which can be 
assembled in a symmetric four tensor g^ui^)- This tensor has a geometrical interpretation 
as the metric of space-time, determining invariant space-time intervals ds in terms of local 
coordinates x^: 

ds^ = g^„{x) dx'^dx". (1.4) 

In the absence of gravity, space-time obeys the rules of Minkowski geometry well-known 
from special relativity. In cartesian coordinates {ct, x, y, z) the metric of Minkowski space 
is 



ri^, = diag(-l,+l,+l,+l). (1.5) 
Then the invariant ds"^ takes the form 

ds'^ = rj^i, dx^dx^ = —c^dt^ + dx^ + dy^ + dz^. (1-6) 

^We use greek letters ^, ^, taking values (0,1,2,3), to denote the space-time components of four 
vectors and tensors. 
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For a test particle, there is a simple physical interpretation of invariant space-time 
intervals measured along its world line, in terms of proper time measured by a clock in the 
local rest frame of the particle. As world lines are time-like, any interval along a world line 
satisifies ds^ < 0, and the proper time interval dr measured by a clock at rest with respect 
to the particle is given by 

c^dT^ = -ds^. (1.7) 
The local matrix inverse of the metric is denoted by g^'^{x): 

g^'{x)g,y{x) = 5^, (1.8) 

where 5^ is the Kronecker delta, and the Einstein summation convention has been used. 
This convention implies that one should automatically perform a full summation over any 
repeated upper and lower indices, like u on the left-hand side of this equation. 

From the gradients of the potentials, i.e. of the metric tensor, one constructs a quantity 
known as the connection, with components F^^^'^ given by the expression: 

= ^ g^"" {d^g.u + d,g^^ - d^g^,) . (1.9) 

At this point this expression may simply be considered a definition. It is a useful definition 
because this combination of gradients of the metric occurs very often in equations in 
general relativity. The deeper mathematical reason for this stems from the transformation 
properties of expressions for physical quantities under local coordinate transformations. 
However, it is not necessary to discuss these aspects at this point. 

A physical argument for the importance of the connection is, that it gives directly the 
gravitational acceleration of a test particle in a given coordinate system by the following 
prescription, due to Einstein: the world lines of point-like test particles in general relativity 
are time-like geodesies of space-time, the latter being considered as a continuous pseudo- 
euclidean spaced with metric g^j,^. Time-like geodesies are lines along which the total proper 
time between two fixed events / dr is an extremum; such world lines are solutions of the 
equation of motion 

^ + r,/^^^o, (1,10) 

A geometrical interpretation of the connection is therefore that it defines geodesies. As 
concerns its dynamical interpretation, there is some similarity of eq. ( |1.1CI| ) with the Lorentz 



force in electrodynamics in that the proper acceleration depends on the four velocity (we 
use an overdot to denote a derivative w.r.t. proper time). However, whereas the Lorentz 
force is linear in the four velocity, the gravitational acceleration depends quadratically 
on the four velocity. The coefficients of this quadratic term are given by the connection. 



^ The term pseudo-euclidean here refers to the fact that the metric is not positive-definite, as it would 
be in a true eucHdean space; space-time has a lorentzian signature (— , -f , -t-, -f ), which imphes that a 
general metric g^n, has one negative (time-like) and three positive (space-like) eigenvalues. 
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which therefore directly determines the force. Notice, that the gravitational acceleration 
does not involve the mass of the test particle. This is the mathematical content of the 
famous equivalence principle, which states that in a fixed gravitational field all bodies fall 
with the same acceleration, independent of their mass. 



Eqs.( l.lO ) and ( |1.9| ) express the force on a test particle in terms of a field (the connec- 
tion), and the field in terms of potentials (the metric). This provides the basis for the field 
theory of gravity formulated by Einstein in his theory of general relativity. 

What is still lacking here of course is a set of equations determining the fields or 
potentials in terms of given sources, like Maxwell's equations allow one to compute the 
electromagnetic fields in terms of charges and currents. General relativity provides such 
equations; they are known as the Einstein equations, and they determine the gravitational 
fields arising from given distributions of mass-energy and momentum. These distributions, 
the sources of the gravitational fields, are described by a symmetric tensor T^jj{x) with 10 
independent components, equal to the number of independent components of the metric; 
T^y is called the energy-momentum tensor. In order to give a concise formulation of 
Einstein's equations, it is useful to first introduce one more geometrical concept, that of 
curvature. 



1.4 Curvature 

Curvature is a important concept, because it characterizes the space-time geometry in a 
way which is basically coordinate independent. In order to illustrate the idea, we first 
present a heuristic discussion. Then we turn to mathematical definitions, which allow 
us to write Einsteins equations in a compact way with a direct and elegant geometrical 
interpretation. 

Curvature is defined in terms of the properties of surfaces; it is a measure of their 
deviation from a flat plane. Curvature can be expressed conveniently in terms of the change 
of a vector when transported parallel to itself around a loop enclosing a two-dimensional 
surface element. In a plane, displacement of a vector parallel to itself will always result 
in the same vector when you return to the starting point: there is absolute parallelism, 
but no curvature. In a curved surface this is no longer true. Consider the example of a 
sphere, a surface of constant curvature. If a vector is moved southward from the north pole 
parallel to itself, along a great circle, it ends up on the equator cutting it at a right angle. 
Now transport it along the equator over an arc of angular extension 6, always keeping at a 
right angle. Finally, move it back to the north pole along a great circle, keeping it parallel 
to the great circle and to itself. It then ends up at the north pole rotated by an angle 9 
with respect to its original orientation. Thus parallel transport of a vector along a closed 
loop on a curved surface does not leave the vector invariant in general. 

In the example of the sphere, the rotation of the vector is proportional to the distance 
traveled along the equator, as measured by the angle 6. This is also proportional to the 
area enclosed by the loop consisting of the two sections of great circles through the north 
pole and the arc of the equator. For a sphere with radius R, the change dV in a vector V 



6 



when the loop encloses an area dA = R'^d6 is 

dV = Vde = ^dA = KVdA. (1.11) 

The constant of proportionality K = 1/ R^ between the relative change in the vector dV/V 
and the change in the area dA is a measure of the intrinsic curvature of the surface, known 
as the Gaussian curvature. Note, that the Gaussian curvature is large when the radius 
of the sphere is small, and that it vanishes in the limit i? — > oo, as one would expect 
intuitively. 

This idea can now be generalized to situations in higher-dimensional spaces. Of course, 
in a space of three or more dimensions, there are many surfaces one can draw through a 
given point. Each of these surfaces may have its own curvature. Our first conclusion is 
therefore, that in general the curvature has many components. 

Fortunately, it is only necessary to consider curvature components in surfaces which are 
linearly independent. To this end, let us choose a locally non-degenerate coordinate system 

and consider in particular the surfaces defined by keeping all coordinates constant except 
two, say x'^ and . 

Let dT^^" = —d'E'^^ be the area of an oriented surface element in the x^-x^-plane which 
contains the point with coordinates (xq). Any vector with components transported 
parallel to itself around the loop defined by the circumference of this surface element now 
changes to first order by an amount 

dV' = -^R^J{xo)V'^di:'^\ (1.12) 
The constants of proportionality R '^{xq) define a tensor known as the Riemann curvature 



tensor at the location (xq). A comparison with eg. ( |1 . 1 1| ) shows, that the Riemann tensor 



is a higher- dimensional generalization of the Gaussian curvature. From eq. (|1.12|) it follows 



immediately, that a space-time is completely flat if and only if all components of the 
Riemann curvature tensor vanish. 

Defining parallel transport of a vector along a line element dx'^ in terms of the Riemann- 
Christoffel connection ( |1.9| ): 

dv^" = -dx^r^fv, (1.13) 

a direct computation for an infinitesimal closed loop gives the result that the components 
of the Riemann tensor can be expressed as 

R,u.^ = d,Tj - d^Tj - [T,,T,l\ (1.14) 

where the square brackets denote an ordinary matrix commutator of two connection com- 
ponents, considered as 4 x 4 matrices (F^^)^'^, the first one with a = fi, the second one with 
a = u. 
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Using eq.( p,.9| ) for the connection, we can ultimately express the curvature in terms of 
the metric and its first and second derivatives. It may then be checked from this explicit 
expression that it identically satisfies the relation 

Rfjiux' + Rux^!^ + Rxpiu'' = 0. (1.15) 

This relation is known as the Bianchi identity. It is analogous to the homogeneous Maxwell- 
equation, which implies that the electric and magnetic fields can be obtained from a four- 
vector potential A^. In a similar way the Bianchi identity for the curvature tensor implies 
that the curvature tensor can be obtained from a symmetric tensor g^j,^. 

We wish to emphasize here, that the curvature of space-time is not a measure for the 
strength of the gravitational fields or forces; rather, it is a measure for the variation of the 
gravitational fields in space and time, the gradients. 



1.5 The Einstein equations 

The Einstein equations are second order non-linear partial differential equations for the 
metric of space-time in terms of a given energy-momentum distribution. As a result, these 
equations describe variations of the fields in space and time, rather than the fields them- 
selves. In view of the role of curvature as a measure for such variations and the geometrical 
interpetation of gravitational forces, it is quite natural that the Einstein equations should 
admit an interpretation as equations for the curvature of space-time. However, rather than 
the full Riemann tensor, it turns out that these equations involve only its contracted forms, 
the Ricci tensor R^^, and the curvature scalar R: 

R>^u = R^xu\ R = R; = g'^" R^.. (1.16) 

Note that we adhere to the general convention that contraction with an inverse metric 
can be used to raise an index of a vector or tensor, turning a covariant component into a 
contravariant one; similarly contraction of an upper index on a vector or tensor with the 
metric itself is used to lower it to make a covariant component out of a contravariant one. 

The reason that only the contracted forms of the Riemann curvature appear in the 
Einstein equations is, that in four-dimensional space-time the Riemann curvature tensor 
has too many components, twice as many as the metric. Therefore fixing the metric 
components in terms of the full Riemann tensor would generally lead to an overdetermined 
system of equations. In contrast, the Ricci tensor is symmetric: R^i, = Ry^, and has 10 
independent components, precisely equal to the number of metric components to be solved 
for. It is therefore to be expected that the Einstein equations take the form of an equation 
for the Ricci tensor. Indeed, Einstein's equations for the gravitational fields can be written 
in the simple form 

Riiv - -^QiJiuR = ^^Mv, (1-17) 

or equivalently: 
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= (T,,-^g,,T), (1.18) 



where T = T^^ = g^^'^T^^. Notice the appearance of Newton's constant on the right-hand 
side of the equations. Of course, in empty space-time T^^ = 0; the second version of 
Einstein's equations ( [L.18| ) then imphes that the Ricci tensor vanishes: 



R^.u = 0. (1.19) 

Space-time geometries which satisfy these conditions are some-times called Ricci-flat. The 
most interesting aspect of these equations is, that an empty space-time can still support 
non-trivial gravitational fields: the Riemann curvature tensor can have non-zero compo- 
nents even if its contracted form, the Ricci tensor, is zero everywhere. 

This observation is important, because it lies at the heart of the theoretical arguments 
for the existence of gravitational waves. Indeed, as we discuss in more detail later on, one 
can find non-trivial solutions of the Einstein equations representing lumps in the gravita- 
tional field which travel at the speed of light and transport a finite amount of energy and 
momentum per unit of volume from one fiat part of empty space to another fiat part of 
empty space. 

There is a modified version of Einstein's equations in which the gravitational field 
itself creates a non-vanishing energy-momentum tensor, even in the absence of matter. 
This gravitational field behaves like a homogeneous continuum with constant density and 
pressure (one of them negative), satisfying the equations 

R,u = -^g,uA, (1.20) 

with A a constant. Such an equation arises if the energy- momentum tensor is 

Tf,, = g^uK. (1.21) 

The constant A is known as the cosmological constant. Contraction of eq. (|1.20|) with the 
inverse metric gives 

R^?^A. (1.22) 

This shows that in the presence of a cosmological constant empty space-time has a constant 
non-zero scalar curvature. Depending on the sign of the cosmological constant, eq.( p..20| ) 
then admits solutions of Einstein's equations corresponding to empty space-times in contin- 
uous expansion or contraction. This modified version of Einstein's equations is especially 
interesting with regard to cosmological applications. 

Of course in the general case there is a contribution to the energy-momentum tensor 
from all matter and energy which happens to be present in the part of the universe we 
wish to describe. The precise form of the energy-momentum tensor depends on the kind 
of matter or radiation which constitutes the source of the gravitational field. However, 
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it is very often justified and convenient to treat matter and radiation as an ideal fluid, 
the energy-momentum tensor of which depends only on the mass- and pressure-density 
and its flow, as decribed by the local four velocity of an element of fluid at location with 
coordinates (x). Denoting the mass density by p(x), the pressure density by p{x), and the 
local four velocity by the energy-momentum tensor then takes the form 



Tf,u = pgf,u + ij> + p) Uf,u^. 



;i.23) 



Comparison with eq.( |1.21| ) shows, that the cosmological constant can be reobtained from 
an ideal fluid with p = —p = A. 



1.6 The action principle 

The Einstein equations ( |1.17| ) can be obtained from a variational principle; that is, they 
follow by requiring an action S to be stationary under variations 6gf^i, of the metric: SS = 0. 
The action for the gravitational field equations, including a coupling to material energy- 
momentum, is 

s = -i_/A^(_£L« + ,..e,„). (1.24) 

Here Qfj,u{x) represents the specific source term, but generally it is not the energy- momentum 
tensor itself; rather, the energy-momentum tensor is given in terms of 0^jy and its covariant 
trace 6 = QJ^ by 

t;. = e^,-^(7^,e. (1.25) 

However, we consider to be the independent variable, hence by fiat there is no contri- 



bution from this tensor under variations with respect to the metric. Moreover, in eq.( p..24| ) 
g{x) = detg^u{x), and the minus sign is necessary because of the lorentzian signature of 
the metric, which implies the product of its eigenvalues to be negative everywhere. The 
proof of the action principle can be stated in a few lines, but it requires some intermediate 
results which we discuss first. 

To begin with, in stead of considering variations of the metric Qf^^, one can equivalently 
use the variations of the inverse metric g'^^. Indeed, because g^xg^" = S^, it follows that 

Sg.xg'" + gf.xSg'" = 0, (1.26) 

or 

Sg^.u = -g^.xSg^'^g.u. (1.27) 

Thus the variation of the metric can be expressed in terms of the variation of the inverse 
metric and vice versa. 
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Next, we need a rule for computing the variation of the determinant g. This is obtained 
from the formula 



logdetM = TrlogM, 



:i.28) 



which holds for symmetric matrices M, as can easily be checked by diagonalizing M and 
using the rule that the logarithm of the product equals the sum of the logarithms of the 
eigenvalues. From this equation we then derive 
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5{-g) = g^''5g,, = -6g^''g,.. 



;i.29) 



Finally, we observe that the Ricci tensor R^i, is an expression involving only the connection 



r^/. Hence 



(1.30) 

with the variation 5R^i, expressible purely terms of the variation of the connection ^J^. 



6R = dg^'-R^, + g^^dR^,, 
;eri 

The total variation of the action now becomes 



6S 



iJd'xV^g 




:i.31) 



6T ^ 

^ fjLV 



The last term is to be interpreted in the sense that we vary R with respect to the metric 
only via its dependence on the connection. However, it turns out that if we insert the 
expression (|1.9|) for the connection, this whole variation adds up to a total derivative. As 
the action principle applies to variations which vanish on the boundary, this term does not 
contribute to 5S. If we then require the variation of the action to vanish, we indeed obtain 
the Einstein equations: 



SvrG 



R, 



:9iivR 



e 



1 



;i.32) 



The equation with the cosmological constant is obtained by making the replacement 



fJ,U 



— » 2A in the action. 



Chapter 2 



Geodesies 



2.1 Curves and geodesies 

In this chapter we discuss how basic information about the geometry of spaces or space- 
times — generically manifolds — can be obtained from studying their geodesies. The 
geodesies of euchdean differentiable manifolds are smooth curves of stationary proper 
length, meaning that to first order the proper length of an arc does not change under 
small variations of the curve. More generally and more precisely, given a metric on a man- 
ifold (euclidean or lorentzian) the geodesies are those smooth curves for which the proper 
interval between any two fixed points V and V through which the curve passes is both 
well-defined and stationary. The interval among the fixed points can be either a maximum, 
a minimum or a saddle point. 

A direct reason for geodesies to be important is that time-like geodesies represent 
the allowed trajectories of test particles in a fixed background space-time. But in general 
geodesies tell us interesting things about the symmetries and topology of a manifold, which 
often have wider implications for physics. 

First we derive the fundamental equation characterizing geodesies. The proper interval 
As [C-pv] of the arc C-pv' of an arbitrary smooth curve C between the two fixed points 
{V, V) is given by the integral 



This holds for all geodesies of euclidean manifolds, and time- or space-like geodesies 
of lorentzian manifolds; A is a continuous parameter on some interval of the real line, 
parametrizing the curve between the points V and V' in any convenient way, provided 
the coordinates x^{\) of points on the curve are monotonic differentiable functions of this 
parameter: the label of any point must be unique and the tangent vector — dx^/d\ 
must be well-defined. The parameter A is then called an affine parameter. Note that the 
integral does not depend on the particular choice of affine parameter: it is invariant under 
differentiable changes of parametrization A — > A' (A). This freedom to choose for A any 




(2.1) 
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smooth parametrization of the curve is often useful, but requires some care in interpreting 
results explicitly containing A. As any invariant geometrical results with physical inter- 
pretation can not depend on the choice of affine parameter, the best way to proceed is to 
remove any A-dependence from equations before drawing physical conclusions. 

The proper interval A s[Cpp'] is stationary if the integral does not change to first order 
under arbitrary smooth variations of the path while keeping the end points of the integral 
fixed: 

6 j d\w[x{\)] = 0, (2.2) 



where the integrand is w[x(A)] = ^J\gf^ui^x'^\- Now after a partial integration with van- 
ishing boundary terms (the end points remain fixed) the variation of the integrand can be 
written 



6w 




Q^u dx 
w dX 



1 dgxjydx^ dx" )^ 
2w dxt" dX dX r 



provided w ^ 0. The sign on the r.h.s. is positive in euclidean spaces or for space-like 
curves in lorentzian space-times, and negative for time-like curves in lorentzian space- 
times. Light-like curves in lorentzian space-time need special consideration, as ds"^ = is 
equivalent to wdX = 0. 

The variational principle ( |2.2D now implies that for any smooth local variation Sx'^ 
the induced variation 6w must vanish if the curve x^{X) is to describe a geodesic. For 
w 7^ we can divide the right-hand side of eq.(^.3D by w once more, and use the fact that 
wdX = ds, the infinitesimal proper interval which is measured in the laboratory frame and 
independent of the choice of affine parameter A. We can also replace the partial derivative 
of the metric in the second term by the connection; recall, that the connection was defined 
in eq.( p..9|) as 

2 ^ \dx^ dx^ dx- )' ^ ' ' 

This expression is actually the solution of the equation 

Dx^ ~ dx^ ~ ^ ' 

stating that the metric is covariantly constant, a condition known as the metric postulate. 
From the equation above it follows that 

^ = r,,«,.. + r,/,.,. (2.6) 

The vanishing of 5w is now seen to imply a differential equation of the form (1.10) for the 
coordinates of a point moving on the curve as a function of the proper interval parameter 
s: 
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-— + VJ — — = 0. 2.7 

ds^ ds ds 

Note, that the above procedure ehminating the affine parameter A is equivalent to a choice 
of rate of flow of the affine parameter d\ equal to the proper interval ds measured in the 
coordinate frame {x^}. If the signature of the manifold is lorentzian and the curve time- 
like this affine parameter is the proper time r, and eq. (|2.7| ) is the equation of motion of a 
test particle of arbitrary non-vanishing mass, in a fixed space-time with affine connection 
V ^J^{x). Note also that the choice of d\ = ds implies an identity 



u 



where again the plus sign holds for euclidean or space-like lorentzian curves, and the minus 
sign for time-like curves in lorentzian space-times. This is a constraint on the solutions of 



the geodesic equation (|2l^ ). 



That this constraint is consistent with the geodesic equation follows by observing that 

is a constant of motion, or more generally a constant of geodesic flow: 

dH dxf" fd'^x'' ^ dx^dx^\ 

Finally light-like curves (also called null-curves) by deflnition obey 

In this case the affine parameter A can not be identified with proper time or distance. As 
eq.(2.11) is to hold on all points of the light-like curve, the quantity 



must be a constant of light-like motion. Therefore eq.( p.lO|) still holds with s replaced by 



A. From this it follows that it is correct to interpret the solutions of the geodesic equation 

+ r./-7r-7r = o, (2.13) 



c/A2 d\ dX 

satisfying the constraint H{X) = as light-like geodesies. 

Remark. Although we have shown that the geodesic equation implies the constancy of 
H{X) for all types of geodesies, the inverse statement is not true. In particular, for non- 
null geodesies (euclidean, and space-like or time-like lorentzian) the more general condition 
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^1 rplX rl'T^ ^rv^^ flT^ 

with B^y = —By^ an anti-symmetric tensor, is sufficient for H to be conserved. This 
happens for example for charged particles in an electro-magnetic field with field strength 

^ IJ,V 

^M^(^) = (2-15) 

fit 

where q/m is the charge-to-mass ratio of the particle. Similar extensions exists for particles 
with spin or non-abelian charges. 

For light-hke geodesies the equation can be further generahzed to read 

1^ + dA = Witt + -^dA • 

with K an arbitrary constant. It should be noted that if k 7^ 0, then H is only conserved 
by virtue of the fact that it vanishes. 

Summary 

The results of this discussion an be summarized as follows. All geodesic curves, in euclidean 
or lorentzian manifolds, satisfy the equation 

dA^ + dA dA = °- P-"' 
where A is an affine parameter. The quantity 

is conserved by the geodesic flow. In physical applications three cases are to be distin- 
guished: 



1. For euclidean manifolds and space- like geodesies of lorentzian manifolds one can take 
d\ = ds, the proper distance; these geodesies are then characterized by the condition 
2H{s) = 1. 

2. Time-like geodesies of lorentzian manifolds can be parametrized by the proper time: 
dX — dr; such geodesies satisfy 2H(t) — —1. 

3. Light-like (null) geodesies admit no such identification of A with invariant measures 
of the arc length in an inertial frame ('laboratory frame'); nuU-geodesics are charac- 
terized by H{X) = 0. 
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2.2 Canonical formulation 



In D-dimensional space-time eq.( |2.17|) represents a set of D second order differential equa- 
tions in tlie affine parameter A providing a covariant description of geodesies. It is possible 
to replace this set by a set of 2D first order differential equations by going over to the 
canonical formulation of the dynamics of test particles, in which the number of indepen- 
dent variables is doubled. Introducing the momentum variables 



the geodesic equation ( |2.17] ) becomes 



dX 



(2.19) 



(2.20) 



The equations ( |2.19| ) and ( |2.20| ) constitute a pair of first-order differential equations equiv- 
alent to the second order geodesic equation (|2.17|) . A powerful result is now obtained 
by observing that these two equations can be derived from the constant of motion H{X), 
eq. (|2.18|) as a pair of canonical hamiltonian equations with H{X) in the role of hamiltonian: 



With this definition 



H{X) 



dxf" dH 



a^'^Pi^Pv 



dpf, 
dX 



dH 



(2.2i; 



f2.22) 



dX dpf. 

Then any function F{x^^, p^) on the 2D-dimensional phase space spanned by the coordinates 
and momenta changes along a geodesic according to 



dF _ dx'' OF dp^ OF _ 

dX dX dx'' dX dp^ ^ ' ^ 



(2.23) 



where the Poisson bracket of two functions (F, G) on the phase space is defined generally 
as 



dF dG dF dG 



(2.24) 



dx'' dpn dpn dx'' 

The anti-symmetry of the Poisson bracket automatically guarantees the conservation of 
the hamiltonian H by 



{H,H} = 0. 

A similar property holds for all quantities which are constant along geodesies: 



dF 

Ix 







^ {F, H} = 0. 



(2.25) 



(2.26) 
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Formal properties of the Poisson bracket include its anti-symmetry and its linearity in each 
of the arguments: 

{F, aid + a2G2} = ai {F, d} + «2 {F, G2} . (2.27) 
Another important property is that it satisfies the Jacobi identity: 

{F, {G, K}} + {G, {K, F}} + {K, {F, G}} = 0. (2.28) 

This property implies, that the Poisson bracket of two constants of geodesic flow is again 
a constant of geodesic flow: if {F, H} = and {G, H} = 0, then 



{{F,G},H} = 0. (2.29) 

As H is itself constant on geodesies, the above results are sufficient to establish that the 
constants of geodesic flow form a Lie algebra, with the Poisson bracket as the Lie bracket. 



The definition of the Poisson bracket (|2.24|) is not manifestly covariant, and therefore 



at first sight it seems to destroy the build-in covariance of the tensor calculus of differential 
geometry, and of general relativity in particular. However, note that in our definition 
the connection ( p.4|) is alway symmetric in its lower indices; therefore we can rewrite the 
Poisson bracket formula as 

{F,G} = V,F—-—V,G, (2.30) 

where for scalar functions F and G on the phase space we define the covariant derivative 

OF 

V,F = d^F + T^Mg^- (2.31) 

This equation preserves manifest covariance for scalar functions F and G; for example, if 
J{x,p) = J'^{x)p^, then 

v,j = {d,r + r^/j") p, = [d^jYp,. (2.32) 

The manifestly covariant form of the Poisson bracket can also be extended to all completely 
symmetric tensors T^^...^^^ of rank n by contracting all n indices with the momentum px to 
obtain the scalar 

T{x,p) = ^T^^^-^^-{x)p,,...p,„. (2.33) 
nl 

Inserting this scalar into the covariant Poisson bracket, the bracket for a particular tensor 
component T^^'"'^" is obtained by taking, in the expression resulting on the right-hand side, 
the coefficient of the n-nomial in the momenta that one is interested in. For anti-symmetric 
index structures a similar result can be achieved using differential forms, or more generally 
Grassmann algebras |ll8| , |l9i pO |. 
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2.3 Action principles 

The canonical formulation of geodesic flow presented here can be cast into the form of a 
standard hamiltonian action principle. Namely, the canonical eqs. (|2.22|) define the critical 
points of the phase-space action 



S= f\\{x -p- H{x,p)) . (2.34) 



By using the hamilton equation 



OH 

dp,,' 



the momentum variable can be eliminated in favour of the velocity of geodesic flow, to give 
the lagrangian form of the action 

S= f d\]-g^,{x)x^xr (2.35) 



We observe, that this is not the action for the geodesies that we started from, eq. (|2.1|) . In 
particular, it is not reparametrization invariant. A comparison shows, that the lagrangian 
in ( p.35| ) is actually the square of w{\) used earlier, eq. ( p.2|) . Clearly, the critical points are 
therefore the same, with the possible exception of light-like curves. Actually, the variational 
principle for the original action was not even well-defined in the case of light-like curves. 

In this section we examine somewhat more closely the relation between the various 
action principles, and in particular the role of reparametrization invariance. It makes clear 
how one can derive the lagrangian and hamiltonian actions presented above within a unified 
framework. 

The key to the general formulation is to incorporate reparametrization invariance in 
an action of the type ( p.35|) |2l| . This is achieved by introduction of a new variable 



e(A), a function of the affine parameter, which acts as a gauge field in one dimension for 
local reparametrizations of the worldline. More precisely, the variable e(A) by definition 
transforms under reparametrizations A ^ A' as 

e(A) e'(A') = e(A) ^. (2.36) 

This is the transformation rule for the square root of a metric in one dimension (on a curve), 
for which reason the variable e(A) is called the einbein. Written in a slightly different 
way, eq.( p.36[ ) states that e{X)d\ is defined to be invariant under reparametrizations. This 



makes it straightforward to write a reparametrization-invariant form of the quadratic action 
(H): 

r.r 1 /"^ ,^ / 1 / .dxf'dx" \ , , 

S[e;x] = dX Ui^-(^)^^ - • (2.37) 
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The last term, proportional to the constant c^, does not depend on the coordinates x'^. 
However, it is reparametrization invariant and its inclusion is important for reasons given 
below. 

Variation of the full action S[e] x] w.r.t. the coordinates x^ reproduces the geodesic 
equation ( ^.171 ) provided everywhere the differentials dX are replaced by edX. This replace- 
ment makes the geodesic equation reparametrization invariant, as edX = e'dX'. On the 
other hand, the substitution e(A) = 1, which breaks reparametrization invariance, makes 
the old and new equations fully identical. Such a substitution is allowed because it only 
serves to define the affine parameter A in terms of quantities measurable by an observer in 
the locally euclidean or lorentzian laboratory frame, like proper length or proper time. Of 
course, such a choice necessarily breaks the reparametrization invariance, but it preserves 
the physical content of the equations. It amounts to a choice of gauge for reparametriza- 
tions in the true sense of the expression, with the result that modulo a constant the action 



S[e; x] reduces to the simple quadratic action (|2.35|) . 

However, the action 5'[e;a;] contains more information, as one can also consider its 
variation w.r.t. the einbein e. Requiring stationarity under this variation gives a constraint 
on the solutions of the geodesic equation, which after multiplication by 2e^ takes the form 

g^^{x)x''x'' + e^c^ = 0. (2.38) 

With e = 1, comparison with the definition of H in eq.( p.9|) shows that this constraint 
expresses the constancy of the hamiltonian along a geodesic: 

2H{X) = -c\ (2.39) 

For euclidean manifolds one can always choose units such that = — 1, but for lorentzian 
space-times = (—1,0,-1-1) corresponds to space-like, light-like or time-like geodesies 
respectively; in the latter case the constant c then represents the speed of light, which is 
unity when expressed in natural units. 

Finally we derive the original reparametrization-invariant geodesic action ( p.l|) from 
S[e; x]. This is achieved by using eq. ( p.38| ) to solve for e, instead of fixing a gauge: 



e = ±,/\^~k^^FJ^\. (2.40) 
Substitution of this expression for e into the action S[e; x] then leads back (modulo a 



possible sign) to the original action (|2.1|). Therefore we conclude that both this action and 



the quadratic action (|2.35|) follow from the same starting point provided by eq.(|2.37|) 



2.4 Symmetries and Killing vectors 

Often space-times of interest possess special symmetries. Such symmetries are useful in 
finding solutions of the geodesic equations, especially if they are continuous symmetries as 
these generate constants of geodesic fiow. This is the content of Noether's theorem applied 
to the motion of test particles and its generalization to space-like curves; see for example 
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refs.||ll||-[p!5[], [ p!8| , |T9|. In this section the role of continuous symmetries is explored in some 
detail. 

A manifold possesses a continuous symmetry, or isometry, if there are infinitesimal 
coordinate transformations leaving the metric invariant. Finite transformations can be 
obtained by reiteration of the infinitesimal transformations, and donot require separate 
discussion at this point. Consider an infinitesimal coordinate transformation characterised 
by local parameters C,^{x): 

x'^ = -i^[x). (2.41) 

For the line element ds^ to be invariant the metric must transform as a tensor; therefore 
by definition 

(2.42) 

On the other hand, invariance of the metric under a coordinate transformation implies 

9',A^)=9,A^). (2.43) 

Combining the two equations ( p.42| ) and ( |2.43| ) to eliminate the transformed metric ^f^^, 
and inserting the infinitesimal transformation ( p.41|) , expansion to first order in leads 
to the covariant condition 

D^^u + D,i^ = 0. (2.44) 

This equation is known as the Killing equation, and the vector fields defined by the solutions 
^^{x) are called Killing vectors. The existence of a Killing vector is equivalent with the 



existence of an infinitesimal coordinate transformation ( |2.41| ) leaving the metric invariant. 

Let us pause for a moment to reflect on the interpretation of this result. If the coordinate 
transformation is considered to be a passive transformation, the new coordinates represent 



a relabeling of the points of the manifold: x^ and x' ^ in eg. ( p. 41 ) refer to the same point 



V. Now consider the point P', which in the new coordinate system has the same value of 
the coordinates as the point V in the old coordinate system: x'^{T") = x^{V). Eq.( p.43| ) 
states that the new metric at the point V has the same functional dependence on the 
new coordinates as the old metric on the old coordinates at the point V. Thus in the 
new coordinate system the neighborhood of V looks identical to the neighborhood of V 
in the old coordinates. In particular it is possible to map the points in the neighborhood 
of V one-to-one to points in the neighborhood of V in such a way that all distances 
(intervals) are preserved. As this is in fact a coordinate- invariant statement, it implies 
a local indistinguishability of the two neighborhoods of the manifold (as long as it is 
not endowed with additional structures). This then is what is implied by a continuous 
symmetry of the manifold generated by the solutions of the Killing equation. 

Now if two neighborhoods of the manifold are similar in the sense described above, 
then this similarity should also extend to the geodesies passing through the points mapped 
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on each other by the symmetry operation. Indeed, as the new metric at the point V' 
has the same functional dependence on the new coordinates there as the old metric at V 
has on the old coordinates, and not just the same value, the connections and curvature 
components there also become identical after the coordinate transformation. The same 
holds for the solutions of the geodesic equations. Thus the symmetry transformation maps 
geodesies through the point V to geodesies through the point V , and the geodesic flow in 
the neighborhood of these points is completely similar. 

The similarity of geodesic structure generated by symmetry operations implies the con- 
servation of certain dynamical quantities under geodesic flow. This is most conveniently 
discussed in the canonical formalism of the geodesic flow presented above in sect. |2.2| . The 
quantity conserved is the generator of the symmetry transformations, obtained by contrac- 
tion of the Killing vector with the canonical momentum: 

J[i] = eix)p,. (2.45) 

This is a scalar quantity, hence its value is coordinate invariant. It generates the infinites- 
imal coordinate transformations through the Poisson bracket ( |2.24| ): 

Sx^" = -^^'{x) = {J[^],x''}. (2.46) 

One can similarly determine the variation of the momentum under the transformation 
generated by J[^]: 

(2.47) 

= d,.m = r/(x)p. 

with T^'^ = d^j,^^ . This represents a coordinate-dependent linear transformation of the 
momentum. 

That the quantity J[^] is conserved now follows from the observation that the hamil- 
tonian is invariant under the phase-space transformations defined in eqs. (|2.46|) , (|2.47| ). 



Indeed, taking into account the canonical equations of motion ( |2.23| ) and ( ^.26[ ), and using 
the Killing equation (|2.44| ) we obtain 



^ = {J[ilH} = 0. (2.48) 

This result was to be expected, as the hamiltonian depends on the coordinates only through 
the metric, which is invariant by definition of an isometry. 

As the Killing equation (|2.44|) is linear it follows that the Killing vectors define a linear 
vector space: any linear combination of two Killing vectors is again a Killing vector. Let 
the dimension of this vector space be r; then any Killing vector can be expanded in terms 
of r linearly independent basis vectors ej(a;), i = 1, ...,r: 



^(ctj) = aiei(x) + ... + ttrCrix). 



(2.49) 
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Equivalently, any infinitesimal symmetry transformation ( |2.41| ) depends linearly on r pa- 
rameters aj, to which correspond an equal number of conserved generators 

J[^] = «! Ji + ... + arJr, where Ji{x,p) = ef (x)p^. (2.50) 

We now show that these generators define a Lie algebra through their Poisson brackets. 
The key observation is, that the Poisson bracket of any two constants of geodesic flow is 
itself a constant of geodesic flow. This is implied by the Jacobi identity for the Poisson 
bracket, which can be written as: 



{{F, G},H} = {{F, H},G}- {{G, H},F}. (2.51) 

Thus, if the Poisson brackets of F and G with the hamiltonian H vanish, then the left-hand 
side vanishes as well, and {F, G} is itself a constant of geodesic flow. 

A second important observation is, that for two constants of motion which are linear 
in the momentum, like the Jj, the Poisson bracket is also linear in the momentum: 

{ Ji, Jj} = (e^jd.et - e'^d.e';) p^. (2.52) 

As the {Ji}l^i are supposed to form a complete set, the expression on the right-hand side 
can be expanded in terms of the basis elements, with the result that with a certain set of 
coefficients f^j^ = —fji^ we have 

{■J^,JJ} = fi.'Jk. (2.53) 
An equivalent statement is that the basis vectors of the Killing space satisfy 

ep^ef-e^d^e^ = /,/ (2.54) 

Eq.( |2.53D , combined with the observations that the Jj form a linear vector space and that 
the Poisson bracket is bilinear and anti-symmetric in its arguments complete the proof 
that the generators of the symmetry transformations define a Lie algebra with structure 
constants f^^j'^. 

Of course the linear p.d.e. (|275^ ) encodes the same information. It can be interpreted 



in terms of the Lie-derivative, which is the operator comparing the value of a function at 
the point V with that at the point V which has the same value of its coordinates after an 



arbitrary transformation of the form (p.41|) as the point V had before the transformation: 



x'^(V') = x^(V). For a scalar 0(x), for which 4>'{x') = 0(x), this gives 

[0] (x) = ling [<f)'{x) - <j){x)] = CdM^), (2.55) 

where the ling[Q] denotes the part of the expression Q linear in ^. For a vector v^{x), for 
which 

v"^ix') = v'^ix)^, (2.56) 
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one finds similarly 

(£5 [v] f (x) = lin^ - vf'ix)] = CduV^ix) - d^^^'v''{x). (2.57) 

For higher-rank tensors the construction of the Lie-derivative gives the result 



(£g [r])^^-^'= (x) = lin^[T'^i-^fc(a;) -T^i-^^x)] 

= ^'^(9^T^i-^*(x) - 9^^^iT^-'^fc(x) - ... - d^i^'^Tt'^-''{x). 



(2.58) 



Now comparing eqs.( p.54"|) and ( p.57| ) we see that the Lie-algebra property of the Killing 

vectors may be expressed as 

^e,N = -C,^[e,] = ^.'^efc. (2.59) 

Finally, in the phase space where by definition d^Pfi = 0, we can also write the transfor- 
mation ( 2.47|) for the momentum as 



Thus we have established a close relationship between Poisson brackets and Lie derivatives, 
and between symmetries. Killing vectors and Lie algebras. 



2.5 Phase-space symmetries and conservation laws 

In the previous section it was established that continuous symmetries of manifolds are 
generated by Killing vectors and imply conservation laws for certain dynamical quantities 
(the generators J) along geodesies. However, not every constant of geodesic flow is neces- 
sarily connected with a Killing vector. This is most easily seen from the example of the 
hamiltonian itself, which is conserved but not of the form (|2.45|) . Rather the hamiltonian 
is quadratic in the momenta. This leads us to consider the possibility of conservation laws 
for quantities which are higher order expressions in the momenta. 

Let J{x,p) be any function on the phase space of coordinates x^(A) and momenta 
Pi_i{\), associated with a geodesic, which is non-singular in the momenta. In others words, 
we suppose it has an expansion of the type 

Here the coefficients J^''^ can be taken to be completely symmetric in all k indices. We ask 
under what conditions such a function can be a constant of geodesic flow: 



(2.62) 
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To obtain the answer, insert the expansion ( p.61| ) into this equation and compare terms at 
equal order of powers of momentum. At zeroth order one finds 

9^/°) = 0, (2.63) 

i.e. J(°) is a constant. Such a constant can always be absorbed in the definition of J and 
is of no consequence. From now on it will be dropped. 

At first order one finds back the Killing equation ( p.44| ): 



These terms we have already discussed in detail. 

Finally, for /c > 2 we find a generalization of the Killing equation for the co-efficients 
j{k) ^lx...^^k(^rJf.y_ 

The parentheses denote full symmetrization over all component indices enclosed. The so- 
lutions of these generalised Killing equations are called Killing tensors. It is now obvious 
that the hamiltonian is included, because the metric g^u is covariantly constant, and there- 
fore automatically a solution of eq. ( p.64| ) for k = 2. Moreover, this solution always exists, 
independent of the particular manifold considered. In addition, for specific manifolds (spe- 
cific metrics (/^^y) other solutions may exist as well. For example, in spherically symmetric 
spaces the square of the angular momentum: 



^ (xY-ix-pf), (2.65) 

is conserved and quadratic in the momenta by construction. Other examples will be en- 
countered later on. 

Just as the Killing vectors, the Killing tensors generate symmetry transformations, not 
of the manifold but rather of the corresponding phase space spanned by the coordinates 
and the momenta (tangent vectors). In particular: 



Sx'^ = {J(x,p),xn = -J«^(x) - Y: 7T-^,J^^^'''-''{x)p,,...p,,. (2.66) 

k=2 y'^ " 

We observe that the transformations generated by the higher-order Killing tensors are 
velocity dependent |18 . 



Finally we address the algebraic properties of the constants of geodesic flow. We have 
already remarked, that the Poisson brackets of two such constants yields another constant 
by the Jacobi identity (|2.51|) . In general the rank of quantities on the left-hand side are 
additive minus one: 

{Z'^),/')} ~ (2.67) 
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Now for / = 1 the rank of jC'^+'-i) is the same as that of J'-'^-*. It follows, that 

(i) for the constants of geodesic flow linear in the momentum, corresponding to Killing 
vectors {k = I = 1), the algebra is closed in itself: it is a Lie-algebra G; 

(ii) for any fixed value of k, the brackets of a J^^^ and a J*^^-* close on the J*^^-*; however, the 
brackets of two elements J^''^ generates an element of order {2k — 1). Hence in general the 
algebra ( p.67| ) is not closed, unless one includes elements of arbitrary rank, i.e. all possible 
k. This implies that once higher-order invariants appear, one may expect an infinite series 
of them. Sometimes this happens indeed, as with the Virasoro algebra for the conformal 
transformation in two dimensions [^]. However, in many cases the series stops because 



the higher order invariants are Casimir-type of operators, or they become trivial in that 
they are just products of commuting lower-order elements (commuting here is meant in 
the sense of vanishing Poisson brackets). 

The Lie algebra G of Killing vectors is characterized by its structure constants /j^- 
As the Poisson brackets of J^''^ with all elements of the Lie algebra spanned by the J*^^-* 
are again of order k, one may choose a basis for the Killing tensors of order 

k, labeled by the index A, and decompose the right-hand side of the Poisson bracket again 
in terms of the corresponding conserved charges jjf'': 

{j1'\J'a'}=9.a'J'b'. (2.68) 

This defines a set of higher-order structure constants. The Jacobi identities then imply, 
that the matrices (Tj)^^ = —Qi^^ define a representation of the Lie algebra ( |2.53| ): 



9,A''9^c'' - 9.a''9,c'' = 4 ' 9k/ ^ nA = 4 ' Tj. (2.69) 

Therefore the constants of geodesic flow 

J^aHx,p) = ^ej)^-^'=(x)p,,...p,„ (2.70) 

span a representation space of the Lie algebra G for the representation {Tj}. The discussion 
of symmetries shows that group theory and Lie algebras can be important tools in the 
analysis of geodesic motion and the structure of manifolds. 



2.6 Example: the rigid rotor 

The above concepts and procedures can be illustrated by the simple example of a rigid rotor, 
which has a physical interpretation as a model for the low-energy behaviour of diatomic 
molecules. As such the example also serves to emphasize the usefulness of geometric 
methods in physics outside the context of general relativity. 

Consider two particles of mass mi and m2, interacting through a central potential V{r) 
depending only on the relative distance r = 1^2 — ril. The lagrangian for this system is 
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L^^rl + r^rl-Vir). (2,71) 

Defining as usual the total mass M = mi + 1712 and the reduced mass n = 17111712/ M, the 
lagrangian is separable in the center of mass coordinates defined by 



MR = miVi + 7712^2, r = r2 — ri. (2.72) 

In terms of {R, r) the expression ( |2.71|) becomes 

L = fR' + ^r^-V{r). (2.73) 

As the center of mass moves like a free particle of mass M, it is convenient to work in 
the rest frame of the center of mass: R = 0. Then we are only left with the problem of 
describing the relative motion of the masses, described by the three coordinates represented 
by the vector r. 

Now suppose that the potential V{r) has a minimum for some separation r = r^, and 
rises steeply for all values of r near this minimum. Then the first excited vibration state of 
the molecule has an energy well above the ground state, and at low temperatures the only 
degrees of freedom that play a role in the dynamics are rotations of the molecule at fixed 
interatomic distance tq. Thus the distance is frozen out as one of the dynamical degrees of 
freedom, and we are effectively left with only two coordinates: the angles {6, (f) describing 
the relative orientation of the two masses; the potential V{r) = V{ro) is constant and may 
be ignored. This leaves as the effective action of the system 

Leff = ^rl (e^ + sinH^^). (2.74) 

This is the lagrangian of a particle of mass /i moving on a spherical surface of radius 
tq. In mathematical physicists' parlance the effective low energy degrees of freedom {9, (p) 
are sometimes called the modular parameters or moduli of the system, and the sphere 
representing the effective low-energy configuration space is called the moduli space. 

Up to an overall scale factor /irg, the lagrangian Lg// represents an action of the form 
( [^.35| ) in a 2-dimensional space with euclidean signature, described by the coordinates 

= 9, x"^ = if and metric 

*=(jsi„%)- P.75) 

As the geodesies on a sphere are great circles, the great circles provide the general solution 
to the equations of motion. The spherical symmetry implies that the system is invariant 
under 3-dimensional rotations (rotations about each of three independent body axes of the 
sphere) and hence there are three Killing vectors generating these rotations. The corre- 
sponding constants of motion are the 3 components of angular momentum L = {Li, L2, L^), 
given explicitly by 
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Li — —sincppe — cot 6 cos (p p^p, 

L2 — cosippe — cot 9 simp p^, (2.76) 
L3 = Pv- 

They span the Lie-algebra so(3) of the 3-dimensional rotation group: 

{Li,Lj} = SijkLk. (2.77) 

As the Poisson brackets of L with the hamihonian vanish, the hamiltonian can depend only 
on the Casimir constant of the rotation group, i.e. the total angular momentum squared L^. 
A direct computation shows, that up to the constant potential the effective hamiltonian is 
in fact proportional to the Casimir invariant: 



(2.78) 



Chapter 3 



Dynamics of space-time 



3.1 Classical solutions of the gravitational field equa- 
tions 



The Einstein equations ( |1.17| ) describe the gravitational fields generated by a specific dis- 
tribution of material sources, as described by the energy-momentum tensor T^,^. In two or 
three dimensions, these fields carry at best a finite number of physical degrees of freedom. 



related to the sources and/or the topology p3[. In four or higher-dimensional space-time 



the gravitational field is a fully dynamical system with infinitely many degrees of freedom. 
In particular, dynamical gravitational fields can exist in the absence of material sources, 
in empty space-time. In the geometrical interpretation of general relativity such fields 
represent the intrinsic dynamics of space-time itself. 

In eq.( |1.12| ) the curvature tensor was defined in terms of the parallel transport of a 



vector around a closed curve. A manifold is called fiat if after parallel transport around an 
arbitrary closed curve in the manifold the image of any vector coincides with the original. 
Therefore a necessary and sufficient condition for a manifold to be fiat near a given point 
is that the components of the Riemann tensor R^^^ vanish there. This condition is more 
general than to require that the metric takes the form 

n 

ds^ = ±dxl + ^ dx'^, (3.1) 

i=l 

with the sign depending on the euclidean or lorentzian signature of the manifold. Indeed, in 
any differentiable manifold one can always find a local coordinate system in the neighbor- 
hood of a given point, in which the line element can be diagonalized to this form, but this 
is not generally the case globally. Even in cases in which it can be done globally, eq. ( |3.1| ) is 
not the most general solution of the Einstein equations for fiat empty space-time, because 
it is still allowed to perform an arbitrary local change of co-ordinates in this metric. The 
fiatness-criterion introduced above, amounting to the vanishing of the Riemann- curvature, 
is more useful as it holds in any co-ordinate system and in any topology. 
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The existence of non-trivial solutions of the Einstein equations with vanishing cosmo- 
logical constant in matter-free regions of space-time is possible because vanishing of the 
Ricci-tensor does not imply a vanishing of all components of the Riemann-tensor. A for- 
mal mathematical criterion for distinguishing solutions of the Einstein equations in empty 
space from those with material sources is through the Weyl curvature tensor, defined as 
the completely traceless part of the Riemann tensor: 



1 

{d-2) 



(3.2) 

Here d is the dimensionality of the manifold. If we similarly define W^j,^, as the traceless 
part of the Ricci tensor: 

W^,^R^,--^g^,R, W^, = W,^, W^/ = 0, (3.3) 

then we obtain a complete decomposition of the Riemann tensor in terms of traceless 
components and the Riemann scalar: 



^ (3.4) 

+ ^pTTY) (^'^'^^^^ ~ Qij-XQuk) R- 

This decomposition allows a classification of the solutions of the Einstein equations ac- 
cording to which of the independent components iW^ni^Xi W^^, R) vanish and which don't. 
In particular, metrics for which only the Weyl tensor is non-vanishing are solutions of the 
source-free Einstein equations R^^, = 0. On the other hand, there also exist manifolds for 
which the only non-zero component is the Riemann scalar R, whilst the Weyl tensor and 
traceless Ricci tensor vanish: the n-spheres xl + xl + . . . + x'l = p^, and all two-dimensional 
surfaces. For the n-spheres this follows essentially from the spherical symmetry, which 
implies the absence of any prefered direction, in combination with the permutation sym- 
metries between the components of the Riemann curvature tensor. 

For two-dimensional manifolds, it is always possible to write the Riemann tensor as 

R^ukx = -^g^nu^nxR = 2 {gtingvx — gvKgnx) R- (3.5) 

For the (torsion-free) riemannian and pseudo-riemannian manifolds we consider in higher 
dimensions (c? > 3), the number of independent components of the Riemann tensor Ri^vkX, 
taking account of all its symmetries, is given by: 
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N[Rf,^^x\ = — • (3.6) 

This is counted as follows: the Riemann tensor is ant i- symmetric in each of the index 
pairs [fiu] and Hence R[^„][kx] may be considered as a square matrix of dimension 

n = d{d — 1)/1. Also as an {n x n)-matrix it is symmetric under interchange of these 
index pairs. As such it has n{n + l)/2 independent elements. Finally the cyclic symmetry 
in the first three (lower) indices of the fully covariant components of the Riemann tensor 



expressed by the Bianchi identity (|1.15|) adds d{d—l){d — 2){d — 3)/4:\ algebraic constraints, 
eliminating an equal number of independent components. Combining these results one finds 
precisely the number N[Rfj^^^x] in eq. (|3.6| ) above. 

For d > 3 the number of components of the traceless Ricci tensor and Riemann scalar 

is 

= ^^^^ - 1. N[R] = 1. (3.7) 

(for d = 2 the Riemann curvature tensor has only one independent component, correspond- 
ing to the Riemann scalar). Therefore for > 3 the number of independent components 
of the Weyl tensor is 

N[W,,,x\ = ^d{d+l){d + 2){d-3). (3.8) 

Clearly this vanishes for d = 3, and is positive non-zero only for d > 4. Thus we find 
indeed that non-trivial dynamical solutions to the source-free Einstein equations can exist 
only in space-times of dimensions d > 4. As the minimal case d = 4 happens to be the 
number of dimensions of our (macroscopic) world, we observe that for this case the Weyl 
tensor and the Ricci tensor each have 10 algebraically independent components. Therefore 
half of the components of the curvature are determined by the material sources, whilst the 
other half describes the intrinsic dynamics of space-time. 

In the geometrical interpretation of the gravitational field and its dynamics, Einstein's 
famous principle of the equivalence of accelerations and local gravitational fields can be 
formulated as the hypothesis that space-time is a manifold, in each point of which one 



can establish a locally flat (cartesian) coordinate system of the form {3A). Of course a 
fully geometric interpretation of all of physics, and also the equivalence principle, breaks 
down in those points or regions where material sources are located^. But also solutions 
of the source-free Einstein equations can violate the equivalence principle locally if there 
are singularities — for example points where curvature invariants become infinite and/or 
the number of dimensions changes. Physics at these points is no longer described, at least 
not completely, by the Einstein equations, and space-time is not necessarily a manifold in 
the strict sense. In some cases it can still be a manifold with boundaries, where specific 



^Attempts to explain matter as a manifestation of physics in higher-dimensional space-time, as in 
Kaluza-Klein models or the more sophisticated string theories, aim to redress this situation. 
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boundary conditions liold; but these boundary conditions then represent external data and 
the gravitational field is not determined by general relativity alone. 

Nevertheless singular solutions are important because they are relevant for physics as 
we observe it in the universe: cosmological solutions of the Friedmann-Robinson- Walker 
type with a 'Big Bang '-singularity, and black-hole solutions describing stellar collapse. The 
singularities in these cases are not necessarily physical, as the short- distance physics may be 
regulated by Planck-scale modifications of the theory, which are presently unobservable. 
However that may be, here we simply accept the possibility of singularities and discuss 
examples of both non-singular and singular geometries arising as solutions of the Einstein 
equations. They illustrate many interesting aspects of gravitational physics. 



3.2 Plane fronted waves 

To show that in four dimensions non-trivial gravitational fields can indeed propagate 
through empty space-time, we now consider a class of solutions known in the literature as 
plane-fronted gravitational waves [^, ^. This name derives from the property that there 
is a co-ordinate system in which these field configurations (i) move with the speed of light 
along a straight line; (ii) have a fiat planar wave front, and {Hi) are of finite duration. 
When discussing the physical properties of these solutions it will become clear that these 
wave-like solutions actually behave more like dynamical domain walls. 

We construct a planar wave solution with the direction of motion along the x-axis, and 
the transverse directions defined by the {y, 2;)-plane. Thus we look for solutions of the 
Einstein equations of the form 



/ -1 



V 











\ 



(3.9) 



The corresponding space-time line-element is 



ds^ = -c^dt^ + dx"^ + f{x,t)dy'^ + g^{x,t)dz'^. (3.10) 

We observe the following properties of these metrics: 

- in the (x,t)-plane space-time is fiat (Minkowskian); 

- in the y- and z-directions space-time is generally not fiat, but the value of the metric 
depends only on the distance along the x-axis and on time, not on the values of y or z 
themselves. Hence in the co-ordinate system chosen the components of these gravitational 
potentials are the same everywhere in the {y, 2;)-plane. In particular they do not fall off to 
unity at spatial infinity, showing that the space-time is not asymptotically Minkowskian at 
large space-like distances. Of course, they share this property with plane electromagnetic 
waves, which at a fixed time also have the same amplitude everywhere in the transverse 
directions; 
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- the metric has no inverse if / = or = 0; then the connections are singular. We must 
take special care when this happens. 

With the Ansatz (|3.9|) for the metric one can compute the components of the connection 
using eq.( |1.9| ). Using the notation x'^ = {ct,x,y,z) for /i = (0,1,2,3), the non-vanishing 
components have been assembled in the following list: 
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(3.11) 



The next step in solving the Einstein equations is to compute the components of the 
Riemann curvature tensor using eq.( |1.14|) . The only independent covariant components 



RfiuXn = R^uxgpn '^o^ identically zero are the following ones: 



-R0202 


— ~^ fftt, 




-R0303 
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= -^ggtt, 


-R1212 


f fxxj 




-R1313 
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= - ggtx-, 
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-R2323 


= fgifxgx 


- \ ftgt)- 







(3.12) 



Before discussing the general solution of the Einstein equations for this class of metrics, we 
first use the above results to answer the question which of the metrics singled out by the 
Ansatz (|3.9|) describe flat Minkowski space. If we require all components of the curvature 
tensor ( p.l2|) to vanish independently, the most general solution for f{x,t) and g{x,t) is 



f{x,t) = a + px + ret, 
g{x,t) = b + qx + let, 



(3.13) 



with the additional constraint 



pq = rl. 



(3.14) 
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Thus metrics of the form ( |3.9| ), and corresponding hne-elements ( p.lO| ), describe a flat 
space-time if / and g are hnear functions of x and t, with the additional constraint ( |3.14|) 
between their co-efficients. 



The components of the Ricci tensor R^i, now follow by taking a covariant trace as in 
eg s. ( p.. 16" ). The non- vanishing ones are 



D 9tt 



r, 73 1 / ftx 9tx \ 

C V / 9 



R 



■22 



R 
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j5 fxx 9xx 

f 9 



f ( fxx — fttj + ^ (^fxdx " '^ft9t 



9\9xx- ^9tt] + T - ft9t 



f 



(3.15) 



In the absence of material sources Einstein's equations require the Ricci tensor to van- 
ish. Assuming / and g themselves not to vanish identically, some manipulation of these 
equations then leads to the following five necessary and sufficient conditions: 

fxx ~^ ftt 0, 

1 

dxx Qtt — 0, 

1 _ (3.16) 

fxQx 5- ftdt — 0, 

gfxx + fgxx = 0, 
9 ftx + fgtx = 0. 



The first two equations are free wave equations in two-dimensional (x, t) space. They imply 
that / and g can be decomposed into left- and right-moving waves: 
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f{x,t) = f+{x + ct) + f_{x - ct), 
9ix,t) = g+{x + ct) + g^{x - ct), 



(3.17) 



where (/_(^), and {f+{i])^ 9+iv) functions only of the single lightcone variable 

^ = X — ct, and 1] = x + ct, respectively. Each of these components then satisfies the wave 
equation by itself, and at this stage they can be taken to be completely independent. The 
only ambiguity here is that a possible constant term can be divided in an arbitrary way 
between the left- and right-moving components of the solution. 

If we now substitute this decomposition of / and g in the third equation ( p.l6|) , and 
assume that the first derivatives of the components with respect to m or f (i.e. /4 and 
g'^) are all non-zero, then the only solutions of the Einstein equations correspond to flat 
Minkowski space. This can be seen by dividing the equation by g'_^g'_ to obtain 



^ = = A 

9+ 9- 



(3.18) 



where A on the right-hand side of the last equality is a finite, non-zero constant. This last 
equality holds because the two previous expressions are functions of different variables, 
hence they can be identically equal only if they are constant. It then follows, that 



(3.19) 



Inserting these results into the last two equations ( p. 161) , addition and subtraction gives 
the equivalent conditions 



{^9 + f)9l = 0, 



{^9 - f) 9- 



0. 



(3.20) 



The solution g'J^ = gl = implies f" = f!' = 0. This leads to linear dependence of / and 
g on X and t: they are of the form ( p.l3|) , ( p.l4| ) and correspond to fiat Minkowski space. 

On the other hand, in case one of the second derivatives g!^ is non-zero, one of the co- 
efficients [Xg ± /) must vanish; note that they cannot vanish simultaneously for non-zero 
/ and g. Hence eq.( |3.2CI| ) gives either 



f = >^9 



and 



4 



n = 0, 



(3.21) 



or 



/ 



-Xg 



and 



r = 0. 



(3.22) 



But combining this with the results (|3.19| ), we find again that all second derivatives must 
vanish and that the geometry is that of flat space-time. 

We conclude that in order to get non-trivial solutions, at least one of the first derivatives 
/± fi'i must vanish. The last three equations ( |3.16|) then imply that, after absorbing any 
constants terms in the remaining components, the only solutions are 
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f{x,t) = f+{x + ct), g{x,t) = g+{x + ct), (3.23) 



or 



f{x,t) = f.{x-ci), g{x,t) = g^{x-ct), (3.24) 

with the additional condition 

f" q" 

^ + — = 0. (3.25) 
/ 9 

Ahhough / and g represent waves travehng either to the right or to the left at the velocity 
of light, the non-linear character of the Einstein equations exemplified by the additional 
condition ( p.25| ) do not allow us to make arbitrary superpositions of these solutions. In- 
terestingly, viewed as a two-dimensional theory in (x, t)-space / and g represent chiral 
bosons. Wherever necessary in the following we chose to work with right-moving fields 
f{x,t) = /-(O and g{x,t) = g-{^) for definiteness. 
An example of a simple solution of eq.(|3.25|) is B: 



f{x,t) = cos k{x - ct) = cos k^, ,^ 
g{x,t) = cosh k(x — ct) = coshfc^. 

Note that at ^ = (2n + l)'K/2k the component / has a zero, whilst the solution for g{x,t) 
grows indefinitely with time or distance. However, this does neither imply that the geom- 
etry is singular, nor that the gravitational fields become arbitrarily strong: the apparent 
problems turn out to be coordinate artifacts. This point is discussed in detail in the next 
section. 

Let us consider a wave which reaches a; = at time t = 0, and extends to x = —L, i.e. 
the solution ( |3.26| ) holds for —L < ^ < 0. If < tt/2 the metric components and g^ 



are strictly positive. For ^ > space-time is flat; we can choose a co-ordinate system in 
which f = g = 1. Note that the wave solution ( p. 26 ) assumes these same values at ^ = 0, 



whilst the first derivatives match as well. As a result the solutions are smoothly connected 
and satisfy the Einstein equations even at the boundary ^ = 0. 

For ^ < —L space-time is again fiat, but now we have to use the general solution ( |3.13| ): 



f{x,t) = a + p{x-ct), 

g{x,t) = b + q{x-ct), [^-ct)<-L. (3.27) 

The co-efficients (a, b; p, q) in this domain of {x, t) values are to be determined again by 
requiring the metric to be continuous and differentiable at the boundary ^ = —L between 
the different domains of the solution. This condition is sufficient to guarantee that the 
Einstein equations are well-defined and satisfied everywhere. Imposing these requirements 
we find 
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cos kL = a — pL, k sin kL = p 

(3.28) 

cosh/cL = b — qL, —ksinhkL = q. 
For / and g in the domain {x — ct) < — L we obtain the result 

f{x,t) = cos kL + k{L + x — ct) sin kL, 

(3.29) 

g{x,t) = cosh kL — k {L + X — ct) sinhfcL. 

Conclusion: we have constructed a solution of the Einstein equations in empty space which 
interpolates between two Minkowski half-spaces, connected by a finite region of non-zero 
curvature, where the gravitational fields have non-zero gradients and test particles are 
accelerated with respect to the initial inertial frame. After passage of the wave the test 
particles have acquired a finite, but constant velocity. This shows that they have been 
accelerated, but the acceleration has ceased after t > {x + L) / c. 

The plane-fronted wave constructed here is a well-defined solution of the equation R^^^ = 
everywhere, including the wave fronts x = ct and x = ct — L, as the contributions of 
the metric components gyy = p and gzz = g^ always compensate each other in eq.(|3.2 



However the Riemann curvature itself, which coincides with the Weyl curvature in this 
case, does not vanish everywhere. Indeed, we find from eq.( |3.12| ) 



-R0202 = -R1212 = ±-Ro2i2 = //", 

-R0303 = -R1313 = ±-Ro3i3 = gg", (3.30) 
-R2323 = 0. 

Inserting the explicit solution for / and g for each of the three regions x > ct, —L < 
X — ct < and x < ct — L one finds that the curvature is zero in the two Minkowski 
half-spaces before and after the wave, as expected, but finite and non-zero in the interior 
of the wave. Indeed, the second derivatives /" and g" have a finite discontinuity on the 
wave front of magnitude =f/c^ at x = ct, and — /c^ cos kL, resp. k'^ cosh kL, aX x = ct — L. 
This finite discontinuity then also appears in the curvature components, where it must be 
attributed to a finite discontinuity in the Weyl curvature, and thus can not be attributed 
to a (singular) distribution of sources. But the Riemann- Christoffel connections, which 
depend only on the first derivatives of the metric, are continuous. Then the geodesies and 
the motions of test particles in the space-time described by the wave are continuous, and 
it follows that the jump in the Weyl curvature presents no physical inconsistencies. 
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3.3 Nature of the space-time 

The planar wave constructed above can be characterised as a region of non-vanishing 
curvature sandwiched between two flat space-times and propagating with the speed of 
light — either in the positive or the negative x-direction — such that one of the two flat 
space-time regions grows at the expense of the other. In this section we establish the 
relation between the two flat space-times, before and after the wave. The analysis shows 
that, although they are of course related by a general co-ordinate transformation, they 
are not related by a simple Lorentz transformation: a set of equidistant test particles 
carrying synchronised ideal clocks and thus defining an inertial co-ordinate system before 
the wave are neither equidistant nor synchronised in any inertial co-ordinate system after 
the wave. In this sense the two Minkowski space-times on either side of the region of non- 
zero curvature are inequivalent: in each the stability group of the vacuum, — the Lorentz 
group of rigid transformations leaving the in-product of four-vectors invariant in each of 
the flat regions — is embedded differently in the group of local co-ordinate transformations 
which constitutes the symmetry of the full theory. 

This situation is also encountered in scalar field theories with broken symmetry, when 
there is a boundary separating two domains with different vacuum expectation values of the 
field. One can even show that in some models these two forms of spontaneous symmetry 



breaking — of a scalar field and of the gravitational field — are connected ||2^. Therefore we 
find a new interpretation for the planar gravitational waves: they act as domain boundaries 
separating different, inequivalent flat-space solutions ('vacua') of the Einstein equations. 
These domain boundaries are dynamical: they are not static but move at the speed of light; 
it is shown below that these particular solutions are block waves of constant curvature. 

Regions of space-time enclosed by such dynamical boundaries clearly grow or shrink 
at the speed of light. If we consider our planar wave solution as an approximation to 
a more general domain-wall geometry with a radius of curvature large compared to the 
characteristic scale of the region of space-time we are interested in, then we can use it to 
describe the dynamics of gravitational domain structure, which has interesting implications 
for the physics of the early universe. 

Having sketched this general picture, I now turn to describe the sandwich structure 
of the planar wave geometry more precisely. A convenient procedure is to introduce new 
co-ordinates = (cT, X, Y, Z) related to the original = {ct, x, y, z) by 

cT = ct + ^, 

X = X + , , 

2' (3.31) 

Y = fy, 
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where in terms of the original coordinates 

M^n = y'ff + z'gg'- (3.32) 

This coordinate transformation takes an even simpler form when expressed in the (nor- 
malized) light-cone variables {U,V): 



U = {cT- X) 1^/2 = {ct - x) 7^2 = M, 

V = {cT + X) /V2 = (ct + x + A) /V2 = v + A/V2. 
Indeed the line-element (|3.10|) now becomes 

(is2 ^ -c'^dt'^ + dx"^ + pdy"^ + g'^dz^ 

= -2dUdV -2K{U,Y,Z)dU^ + dY^ + dZ\ 



(3.33) 



(3.34) 



with 



f" n" n" , 



K(U,Y,Z) = - ±-Y^ - ^ = — {¥' - Z'] . (3.35) 
/ 9 9 ^ ' 

Both for t < x/c and t > {x + L)/c we have /" = g" = 0, and therefore the line element 
( ^.34| ) is manifestly flat Minkowskian on both sides of the planar wave. However, inside 



our particular wave solution one has 

f" n" 

■' ^ = - P. (3.36) 



/ 

Then K = k'^{Y'^ — Z'^), independent of U. The coordinate system with the line-element 
( p.34| ) is the closest one to an inertial frame at all times for an observer at rest in the origin 



one can construct. Indeed, the point x = y = z = corresponds toX = Y = Z = 
at all times, and both before and after the wave the X^^ define a global inertial frame. 
Furthermore, inside the wave the new coordinate system deviates from an inertial frame 
only for large values of {Y^ — Z^). We also observe that in this coordinate system the line 
element ( ^.34| ) never becomes singular, showing the previous zero's of detg, corresponding 
to zero's of f{u), to be coordinate artifacts indeed. 

Computing the various geometrical quantities: metric, connections and curvature com- 
ponents, in the new coordinates produces the following summary of results: 

- The only non-vanishing components of the inverse metric are: 

guv^gvu^_^^ gvv = 2K = 2k\Y^-Z^), 

(3.37) 

g^y = g^^ = 1. 

- The only non-vanishing components of the connection are: 
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p Y _ -p V _ -p V _ nl.2Y 



(3.38) 

The only non- vanishing covariant components of the curvature are 



RuYUY — —Ruzuz — — 2fc^. (3.39) 

Thus in these coordinates the curvature components are constant inside the wave, and 
we conclude that our particular planar wave solution is actually like a blockwave: at the 
boundaries f/ = 0, [/ = —L the curvature components jump from zero to the values 
( ^.39| ) and back to zero, whilst in between they are constant. This makes it easy to 



visualize how our particular solution ( |3.26| ) can be used to construct general planar wave 



solutions with polarization along the Y- and Z-axes: essentially it amounts to glueing 
together sufficiently narrow block waves of fixed curvature to obtain solutions in which the 
curvature components vary according to some arbitrary profile. 



The result (|3.39|) also allows one to verify once more by a one-line computation that 



the Ricci tensor vanishes identically for these planar-wave geometries. 

With these results in hand, our next aim is to establish the relation between the flat half- 
spaces before and after the wave. We have already observed that the origin x = y = z = Q 
becomes the origin X = Y = Z = {] for all times t or T, respectively. Next observe. 



that the coordinate transformation ( |3.31| ) is not a Lorentz transformation: the coordinates 
(t, x, z), defining an inertial frame before the passage of the wave, represent a non-inertial 
frame after passage. As a check, note that the jacobian of the transformation is 



= fg^l^ fore = ct-x^O. (3.40) 

Indeed, a system of synchronized, equidistant clocks at rest in the first coordinate system, 
which operationally defines an inertial frame before the wave arrives, is neither synchro- 
nized nor equidistant in the final inertial frame. Moreover, in general it cannot be made 
so by a simple Lorentz transformation: the clocks have been accelerated with respect to 
each other by the wave, and ultimately define a non-inertial frame at late times. This is 
verified by the explicit calculation of geodesies in the next section. 

From this analysis it will be clear, that in principle it is possible to determine the 
passage, and measure the characteristics, of a planar gravitational wave by measuring the 
acceleration it imparts to test masses in the laboratory. Indeed, such measurements can 
determine the functions {A,f,g) from the relation between (t,x,y,z) and (T, X,Y, Z). 
Although in this form the statement holds only for waves of the type discussed, the princi- 
ples for other types of waves are the same, even if the functions characterizing the relation 
between the initial and final Minkowski spaces are generally different. 
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3.4 Scattering of test particles 

In order to study the response of test masses to passing gravitational waves, it is necessary 
to solve for the geodesies of space-time in the gravitational field generated by the wave. 
We are particularly interested in the worldlines of particles initially at rest in a specific 
inertial frame (the laboratory). Observing their motion in the final inertial frame gives the 
scattering data, from which the properties of the gravitational wave can then be inferred 
relatively easily. 

In principle it is very easy to find geodesies: they correspond to the worldlines of 
particles in free fall. Consider therefore a particle at rest in the initial Minkowski frame; 
its four co-ordinates are 



x^{t) = (cr,ei,6,e3), (3.41) 

where the three vector ^ = (^1,^2,^3) fixing the position of the particle is constant; note 
also that by definition for a particle at rest in an inertial frame the co-ordinate time equals 
the proper time r. 

In the co-ordinate system we use, a particle with co-ordinates ( p.41|) is actually always in 



free fall, its co-ordinates in the x'^-frame do not change if it enters a non-trivial gravitational 
field. Of course, in an arbitrary gravitational field the x^-frame is not an inertial frame, 
and even if at late times space-time is flat again the particle will have acquired a finite, 
generally non-vanishing four velocity. However, its position in the final Minkowski frame 
can be found straightforwardly by applying the co-ordinate transformation ( ^.311 ). 



Let us first confirm the above statements by a simple mathematical argument. Then we 
turn to apply the result to find the actual motion of the test particle in the final Minkowski 
frame. We begin by noting that in the x^-frame the metric for the gravitational wave is 
always of the form (|3.9| ), and therefore the only non- vanishing components of the connection 
are those given in eq.( p.ll| ). In particular, in the x'^-frame the components Foo'^ always 



vanish. 

On the other hand, for a particle with the worldline ( pj.41[ ), the four velocity is purely 
time-like and constant: 

= ^ = (c, 0,0,0). (3.42) 
Then the four acceleration in this frame vanishes identically, as does the gravitational force: 

r^/w^u^ = 0. (3.43) 

As a result the geodesic equation ( |1.10| ) is satisfied, proving that ( p.41| ) represents the 
worldline of a particle at all times. 

What we are particularly interested in, is the position and velocity of the particle at the 
moment it comes out of the gravitational wave. It has already been remarked that the 
origin ^ = remains at rest and coincides with the origin of the final Minkowski frame. 
Therefore at late times the origin has laboratory coordinates 
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X^ir) = (cr, 0,0,0). (3.44) 

Let us compare this with the four co-ordinates of a particle away from the origin in the 
(y, 2)-plane: $, = (0,^2,^3) 7^ 0. For such a particle the values of the metric components 
at late times are 

^ l'"""^' OT>L. (3.45) 
5,(0, T) = b - qcT, 

Here {a,b;p,q) have the values quoted before, eq.( p.28| ). For A one finds in the same 
domain ct > L: 

A = ape2 + bqil - (fil + q'il) CT. (3.46) 
As a result, the Minkowski co-odinates of this particle at late proper times are 



Xf(r) = (cTi,Xi,ri,Zi), (3.47) 



with 



cTi = -i {apil + hqil) + f (2 + p\l + q'il) , 

X. = -I {ape, + bq^) + f {pr. + q'^l) , ^3_4g^ 

Zi = b^3 - qizCT. 

Note, that clocks which were originally synchronized — same value of t = r at early proper 
time — are no longer so after passage of the gravitational wave: Ti is not equal to Tq = r at 
late proper time. Before we can calculate the position of the test particle at synchronized 
co-ordinate time T = r, in the final Minkowski frame, we first have to correct for this 
gravitationally induced position-dependent time shift. 

To determine the instantaneous proper distance between the particle with worldine 
Xf (r) and the origin, we have to ask at what proper time f the time co-ordinate Ti equals 
the time Tq measured by a clock at rest in the origin: 

Ti(f) = r, (3.49) 

with the result 

2cr + anfn + bq^'i 

cf = ^ ■ 3.50 

The position at this proper time, = (f), can now be evaluated: 



41 



2 + p2^2 + ^2^2 



V c ( (»g ~ bp) q^l + 2 (g - per) ^ 

- ^ f -jaq- bp) pil + 2{h- qcr) \ 
2+p^e2+q'e3 J' 
Finally we can also determine the velocity of the test particle: 

- = 2.L 2.2 {P'^i + q\l -2^2, -2qi,) , (3.52) 

c 2 + PX2 + q ^ ' 

and result the kinetic energy is 



2 (2 + p^el + 



From this analysis we draw the conclusion, that not only in general test particles are 
accelerated by the gravitational wave, but also that the acceleration and hence the final 
velocity depend on the position of the particle. Therefore no Lorentz transformation with 
any fixed boost velocity v can bring all test particles initially at rest w.r.t. one another 
back to rest simultaneously after passage of the wave. This result explicitly demonstrates 
the inequivalence of the initial and final Minkowski space-times. 



3.5 General planar waves 

Having studied in great detail a very specific planar wave solution of the Einstein equations, 
we now briefiy consider the general planar wave solution defined by metrics of the form 
(PI): 

ds^ = -2dUdV - 2K{U, Y, Z)dU^ + dY^ + dZ^. (3.54) 



It was mentioned below eq. ( p.39[ ) that solutions of this type can be found by glueing 



together block waves of constant curvature. The general solution is obtained analytically 
by allowing for arbitrary dependence of K{U, Y, Z) on its arguments. The expressions for 
the corresponding connection components are 



^uu = ^,u, ^Yu = ^uu = K.Y, Tzu = ^uu = ^,z, (3-55) 



all others being zero. Up to permutation of indices, the only non-vanishing components of 
the Riemann curvature tensor then turn out to be 
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RuYUY — —KyYi RuZUZ — —KzZi 

(3.56) 

RuYUZ = RuZUY = —K^YZ- 

As a result the Ricci tensor takes the simple form 

Rmn = Ruu^M^N^ (3.57) 

with 

RUU = - {Kyy + K^ZZ) = -^transK. (3.58) 

Here ^trans denotes the laplacian in the transverse y-Z-plane. It is convenient to introduce 
complex co-ordinates in this plane: 

C = (r + iZ) /V2, C = (^^ - ^Z) /V2. (3.59) 

The Einstein equations in empty space-time then reduce to the linear second order wave 
equation for K: 

- AtransK = -K^^^ = 0, (3.60) 

The general solution of this equation is a sum of holomorphic and anti-holomorphic func- 
tions of {(, Q. Expanding in Fourier modes it can be written as 

KoiU, C, = E y_ (5„(fc)e-^ V + ^-n(fc)e^'^C) • (3.61) 



n=0 



From eqs.( p.5(^ ) it follows that the constant and linear terms in and C, correspond to 
a vanishing Riemann curvature tensor and hence describe fiat space-time. Therefore it 
should always be possible to find a coordinate system in which these terms are absent. In 
this system the expansion ( p.61[ ) in Fourier modes then starts at n = 2. 



The proper-time equations of motion for a point mass take the form 
dr'^ dT'^ dr dr 



d'Y 2^ d'Z 



(3.62) 



z- 



dr^ ' ' ' dr^ 

Here we have already inserted the conservation law for the proper ^/-velocity 

dU 

— — = 7 = constant. (3.63) 
dr 

The origin of proper time can be chosen on the light-cone ?7(0) = 0, with the result that 
U = jT. It is then possible to exchange the proper-time r for the light-cone time U, which 
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sometimes leads to simplifications. Next, the conservation law for the particle hamiltonian, 
generating time-like geodesic flow, gives 

We can combine the last two equations to find equivalent expressions for the proper-time 
derivatives of (T, X). Adding the square of the X-component of the proper velocity to the 
left- and right-hand side of eq.( p.64| ) gives 





(3.65) 



Defining the velocity in the laboratory frame a.s v = dR/dT, eq.( p.65| ) leads to the result 



dT 1 - 2^^K/c^ 

^ = \\ ■ ^^-^^^ 

This establishes the physical interpretation of K as the gravitational potential determining 
the gravitational contribution to the time-dilation (redshift). We can invert this relation 
to establish an expression for the constant 7 in terms of the gravitational potential K and 
the velocity components: 

c ^ dr ( c ) I 1 — l(? ^ ^ 

An immediate corrolary of this result is that the quantity 

is a constant of motion. This result takes the place of conservation of energy for test 
particles in the field of the gravitational wave. It results from the light-cone co-ordinate V 
being cyclic. 

Thus the equations for (f/, ^) can be integrated in general for arbitrary K(IJ ^Y^Z\ 
Then it only remains to solve the equations for (F, Z) on the second line of eq.( p.62|) . Of 
course, the explicit solutions of these equations depend on the specific form of K(\J, Y, Z). 
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3.6 The gravitational field of a light wave 

So far we have considered gravitational waves in empty space-time. But in nature grav- 
itational waves will most likely be associated with other physical phenomena, with the 
dynamical behaviour of matter and energy. In fact we expect gravitational waves to ac- 
company almost any time- varying distributions of matter and energy, provided the motions 
are not too trivial (like uniform acceleration). 

In this respect is instructive to study a relative simple and well-known process, which 
is nevertheless quite non-trivial: the propagation of electro-magnetic waves. Waves in the 
electro-magnetic field provide an example of a time-varying energy and momentum distri- 
bution in space-time, accompanied by a co-traveling gravitational wave field. The physical 
presence of this gravitational field traveling at the speed of light can be demonstrated 
in principle by considering the motion of an electrically neutral point-like test particle. 
Although it is totally insensitive to the electro-magnetic wave, it does feel the associated 
gravitational wave and will therefore start to oscillate. These orbits of neutral point-masses 
are quite non-trivial; in general they are non-periodic, except for the case of zero-frequency: 
a pair of constant crossed electric and magnetic fields generates pure harmonic motion of 
a neutral point mass by its accompanying gravitational field. This is a quite unexpected 
and strange phenomenon, which may best be compared with the Josephson effect in su- 
perconductivity, in which a constant potential generates an oscillating current. In fact we 
will demonstate an even closer analogy in the case of scalar fields in the next section. 

With the conventional energy-momentum tensor for the electro-magnetic field, the 
Einstein-Maxwell equations describing the coupling between gravity and electro-magnetism 
read |23 



= {f,,F/ - \g,^F' ) , (3.69) 

and 



D^F^^ = 0. (3.70) 

Because the energy-momentum tensor of the electro-magnetic field is traceless, the first 
equation implies that the scalar curvature vanishes, R = 0. 

A free monochromatic plane electromagnetic wave of angular frequency uj = kc and 
traveling in the x-direction is described in the Coulomb gauge by a vector potential 

A = asmk{x — ct)), (3-71) 

where a is a constant vector in the y-z-pla.ne: = 0. The corresponding electric and 
magnetic fields are 



OA ^ ^ ^ ^ 

E = — — = ujd cosk{x — ct), B = VxA = kxd cosk{x — ct). (3.72) 
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The wave-vector k = {k, 0, 0) is orthogonal to a, hence 



k X a 



= \ka\ = ka, and 



E 



cB 



Eocosk^x — ct), Eq = ua, (3.73) 



whilst E ■ B = 0. The wave (|3.71|) remains a solution in general relativity provided the 
covariant derivative in eq.( |3.7U| ) acts as an ordinary derivative. Now from the symmetry 
of the Riemann-Christoffel connection F^^'^ = T^^J^ and the anti-symmetry of the electro- 
magnetic field-strength tensor F'^'^ = —F'^^ it follows that the covariant divergence of 



eq. (|3.7Q|) can be written explicitly as 



D^F^"" = d^F"'' + T^^^F^'' = 0. (3.74) 

Therefore the covariant derivative turns into an ordinary derivative if the contracted con- 
nection vanishes: T^^^ = 0. This happens if the metric satisfies g = detg^i, = constant. 

The symmetry of the plane-wave configuration implies that the energy momentum 
of the electro-magnetic field varies only in the x- and t-directions; thus we expect the 
curvature to change only in those directions as well. Moreover, to keep the electric and 
magnetic components of the wave perpendicular we require the metric in the y-z-p\ane to 
be trivial. Finally, as the plane-wave configuration depends only on the variable {ct — x) 
it is convenient to switch to light-cone coordinates 

ct-x ct + x 

The above requirements can be satisfied by taking a metric of the form 

ds'^ = -2dudv - 2K{u,y,z)du'^ + dy"^ + dz^ . (3.76) 

This Ansatz for the metric is the same as that for free gravitational waves, eqs. (|3.34|) , 
( ^.54| ). It has constant determinant g = —1 and is manifestly flat in the transverse plane. 

The non-vanishing components of the Riemann-Christoffel connection are the same as 
in eq.( p.55D : 



^ UU ~ ^^,U) ^ yu ~ ^ UU ~ ^^,y} ZU ~ UU ~ ,Z ' \^ ' ' ' ) 

From these expressions one can explicitly verify that the contracted form of the connection 
vanishes: T-^^^ = 0; this implies that Maxwell's equations are satisfied indeed. 
As before, the non-zero components of the Riemann curvature are 



Ryuyu ^,'yyi -Rzuzu -^jZzi 



Ryuzu Fizuyu ,yz ' 



(3.78) 



Then the Ricci tensor reads 



R^.u = RuuS;5:, (3.79) 
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where again 

Run = - {K,yy + K,,,) . (3.80) 

From eg. ( ^.721 ), the energy-momentum tensor of the electro-magnetic wave takes the form 
- \g,A = {fL + Fly) 6;5: = 2k'a' cos' {V2ku) 5^5:, (3.81) 

with 

Fui = Ai^u = ^ {Ei + eijcBj) = V2kai cos V^ku, (3.82) 

cv 2 

with eij the 2-dimensional anti-symmetric epsilon symbol in the transverse {y, z)-plane. 

Comparing egs. (|3.81|) and ( p. 791) , and introducing as in eg. (p.59|) the complex coordi- 
nates C = {y + W^) I \/2 and C = {y — iz) / v^, we see that they match indeed provided 
the single condition: 



E' + c^m = El cos' V2ku, (3.83) 



is satisfied. Remarkably, this eguation is linear in the metric component K. Therefore its 
general solution can be given as 



K{u,y,z) = ^^CCEl cos'{V2ku) + f{u;0 + />;C), (3.84) 



where the first term on the right is a special solution of the inhomogeneous Einstein egua- 
tion, whilst f{u] is an arbitrary complex analytic function of u and (, representing the 
general solution of the homogeneous eguation |2^. Thus the general solution of Einstein's 



eguations is seen to be a superposition of a special gravitational wave, generated by the 
energy-momentum of the electro-magnetic field, and an arbitrary free planar gravitational 
wave traveling in the same direction. 

It is straightforward to extend the solution to the case of a pulse of radiation described 
by a vector potential 

A{x,t) = J dk dk sin k{x — ct) , (3.85) 

where the direction of polarization determined by dk may depend on k provided it is 
transverse to the axis of propagation: k ■ dk = 0. The energy momentum tensor of the 
pulse is similar to that of the monochromatic wave (|3.81|) : 

FfixF^ — -^QfiuFf^x — (^F^^ + F^y) S^Sjj 

(3.86) 

= 6^6^ ^ J dkj dqEk-Eg (cos V2{k + q)u + cos V2{k - q)u 
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The corresponding metric is again of the form ( |3.76| ), with 

Aire G ^^'^^^ 

The physical nature of the special gravitational wave generated by the electro-magnetic field 
can be investigated by solving the equations for time-like geodesies, representing the orbits 
of electrically neutral massive test particles. Parametrizing the geodesies by the proper 
time T and using an overdot to denote proper-time derivatives, the geodesic equation is 

xf" + T^/x^x^ = 0. (3.88) 

Consider the simplest case of a monochromatic electro-gravitational wave, with f{u; () = 0; 
then the metric component K can be written as 

K,m{u, x,y) = ^ {l + cos 2V2ku) (y^ + z^) . (3.89) 
In this expression we have defined the constant 



nl = '-^El. (3.90) 



Using this in the connections (|3.77|) , the first geodesic equation for u{t) becomes 



M = 0. (3.91) 
It follows that the momentum ?i is a constant of motion: 

M = 7. (3.92) 

In addition we can write down two other constants of motion. The first one is the hamil- 
tonian 

H = jg.^x^x'^ = (3-93) 

where we have used that for time-like geodesies the geodesic interval is measured by proper 
time: ds'^ = —c^dr"^. From the expression for the metric in eq. ( p.76| ) it then follows that 

V = ^ (c^ + y'^ + z^) - ^K. (3.94) 
27 ^ ' 

The other constant of motion is the angular momentum in the z-direction 



J = yz - zy. 



(3.95) 
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It remains to solve for ?/(r) and z{t). The geodesic equations for these coordinates take 
the form 



y = - ^ (1 + COS 2cjir) y, = - ^ (1 + cos 2cjir) z, (3.96) 



with the frequencies ujq and cJi defined by 



UJo 



7^072, tui = 7A;V2. (3.97) 



In the limit of static electric and magnetic fields, k = ui = with Eq and Bq fixed, they 
reduce to ordinary oscillator equations, with purely harmonic solutions of frequency c^o- 
Hence we find that a neutral test particle in static crossed electric and magnetic fields 
executes a harmonic motion in the y-z-plane, as a result of an induced gravitational wave 
| 28| . As pointed out in the introductory remarks, this phenomenon as well as a similar 



one for the gravitational effect of a scalar field with constant gradient p7| , |28[| , recall the 
Josephson effect in superconductivity. 

For uji ^ eqs. (|3.96|) are particular instances of the Mathieu equation, the solutions of 
which are well-known [^. In general they are not periodic. For the present application 
it is convenient to proceed as follows. Introduce the new affine parameter A = cjit; then 
writing rj = ujq/uji for the ratio of the two frequencies we obtain the equation of motion 

^ + (1 + COS2A) y = 0, (3.98) 

and a similar equation for -2(A). In practice we always have r/ <^ 1 therefore we can expand 
y{X) in powers of r/^: 

oo 

y{X) = Y: Cn{\W\ (3.99) 

n=0 

This can be solved recursively for the c„(A): 

The solution of these equations is straightforward. We fix the constants of integration such 
that ?/(0) = yQ = Co(0) and y'{0) = Py = c'q{0), and therefore 

Cfc(O) = 4(0) = 0, k = 1,2,... (3.101) 

Here the prime is used to denote a derivative with respect to A. The solutions then take 
the form 

yiX) = yoa{X) + (3yb{X), 

(3.102) 

^(A) = zoa{X) + P,b{X), 
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where the universal functions a(A) and 6(A) contain all the dependence on the affine pa- 
rameter A. They possess the expansions 

2 4 

a(A) = 1 + ^ai(A) + ^a,{X) + ... 
4 lo 

(3.103) 

6(A) = A (^1 + ^6i(A) + 1^ 62(A) + 

with 



ai(A) = + I cos2A, 

25 A2 A^ /3 An 1 
a2(A) = + X y + ( 4 " y ) """^ ^ + ^ ^'"" ^ + 32 ' 

6.(A) = ^-y-^ + ^os2A, (3.104) 

45 A2 A^ /5 sin2A /5 A^X 

62(A) = + - + AM — — + cos2A 

' 32 4 30 V2 y 2A V4 6 J 

3 sin4A 1 

^ \ cos4A. 

16 4A 32 

A non-trivial check on this result follows from the conservation of the x-component of 
angular momentum. Its value is fixed by the initial positions and velocities as 

J = VoVz - ZoVy, (3.105) 
where Vi = uJiPi. It is constant because the functions a(A) and 6(A) satisfy the relation 

a 6' - 6 a' = 1. (3.106) 
This may be checked order by order in powers of r/^. In particular, one finds 

ai - Xa[ + 61 + A6; = 0, 

(3.107) 

a2 — \a'2 + 62 + A62 + aibi + A {aib[ — bia[) = 0. 

Thus the conservation of angular momentum relates the functions a(A) and 6(A). 

The expansion ( |3.99| ) assumes rj < 1. From the explicit results for a(A) and 6(A) we 
can infer that the coefficients remain small for values of A such that Xr] = ujqt < 1, i.e. for 
proper times r < I/uq. Numerically 
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uo (rad/sec) = 0.29 x 10"^^ Eq (V/m). (3.108) 

c 

For non-relativistic motion 7v^ «i c, whilst the highest electric fields in the laboratory 
are generated by lasers, with Eo ~ 10^^ V/m [^]. This limits c^o in practice to values 
less than 10^"^ rad/sec, and therefore almost always ujq ^ uji. Note that on top of the 
slow oscillation we observed for the case of static (zero frequency) electro-magnetic fields, 
there is now a modulation with frequencies 2nciJi, where n is an integer number, but unless 
Ui ~ ujq the amplitude of these oscillations is very small. 

A theoretical upper limit on the frequency ujq is provided by quantum electro-dynamics, 
as the vacuum breaks down due to electron-positron pair creation for values of the order 
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E^„, = ^^^^ = 4.1 X lO^^V/m. (3.109) 

eXe 

Here mg is the electron mass, and Ag its reduced Compton wave length; for details we refer 
to ref.|^. The critical field strength corresponds to a frequency upper hmit of cuq ~ 1 
rad/sec. 

From the solution ( p.l02| ) of the equations of motion it follows that a point mass at 



rest in the origin remains at rest at all times. Now clearly, the choice of origin in the 
y-z-pla.ne is arbitrary due to translation invariance. Therefore the non-trivial nature of 
the solution becomes manifest only if we compare the relative motion of test particles. 
Fortunately, using one test particle initially at non-zero [yo, zq) and one at rest in the 
origin, our solutions are seen to determine the relative acceleration between these point 
masses. As a result they provide a true measure for the curvature of space-time generated 
by a light- wave. 



3.7 Symmetry breaking as a source of gravitational 
oscillations 

The planar gravitational wave solutions described earlier in this chapter as block-waves, 
arise as transition regions between fiat domains of space-time, which are related by a 
coordinate transformation that cannot be reduced to a (linear) Lorentz transformation. 
We have therefore interpreted them as domain boundaries between regions of space-time 
trapped in inequivalent classical vacua of the gravitational field (in the absence of a cos- 
mological constant). This is possible because fiat space (zero curvature) is associated with 
a value of the gravitational field, the metric g^^iy, which in the appropriate Minkowskian co- 
ordinate system ('Minkowski gauge') is constant, non- vanishing and Lorentz- invariant, but 
not invariant under general coordinate transformations. In this sense the general coordi- 
nate invariance of the Einstein equations is broken spontaneously by the vacuum solutions 
representing fiat Minkowski space-time. 
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A domain structure arising from the presence of degenerate but inequivalent classical 
vacua is a well-known consequence of spontaneous symmetry breaking in scalar field theo- 
ries. In this section we take up this aspect of the the field theory and show, that dynamical 
domain boundaries in scalar field theories (scalar-type shock waves, if one prefers) are also 
sources of planar gravitational waves of the kind analysed in detail above This is 

another instance of the general observation that almost any type of dynamical situation 
with time-dependent mass-energy distributions is accompanied by gravitational fields. 

The simplest scalar field theory with a continuously degenerate set of classical vacua is 
that of a complex scalar field with a 'mexican hat' type of potential, as described by the 
lagrangian 
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C = -g^''d,ct>*d,<P - 0*0-^^ , (3.110) 



In this lagrangian hc\ is a dimensionless coupling constant and /i is an inverse length 
{fih/c is a mass). In the natural system of units in which ^ = c = 1 we can take A to be 
dimensionless and /x to represent a mass. In the following equations one then can substitute 
c = 1 everywhere. 

This lagrangian ( |3.110|) is invariant under rigid U{1) phase shifts — > e*°0, with a 



constant. As a result for /i^ > there is an infinitely degenerate set of solutions of the 
equations of motion minimizing the potential (refered to as classical vacua), of the form 

0=-^e^ (3.111) 

with 6 an arbitrary real constant. Clearly the solutions themselves are not invariant un- 
der the U{1) transformations: the classical vacuum solutions break the U{1) symmetry 
spontaneously. However, because of the U{1) symmetry of the dynamics these solutions 
are physically indistinguishable, and the absolute value of 6 is unobservable. On the other 
hand, variations in the vacuum angle between different regions of space-time can have 
observable effects [031 EM. 



Parametrizing the complex scalar field as 0(x) = p(x)e^^(^VV2, the la grangian becomes 
1 „„. . X f o 2u^V 1 



^ = -7^g''d,pd,p - - - ^ - -p'g^^d.ed^e. (3.112) 



This shows that in the spontaneously broken mode of the U{1) symmetry the theory 



desribes actually two real scalar fields: the field x = P^/^v ^Z-^' obtained by shifting p by its 



classical vacuum value nJ2/X, with a mass (or inverse coherence length) equal to /i, which 



can become arbitrarily large; and the strictly massless field a = fi6J2/X, representing the 



Goldstone excitations arising from spontaneous symmetry breaking ^ . 

Clearly at energies small compared to the scale fi set by the mass of the field x the 
excitations of this massive field play no role and one only has to consider the Goldstone 
field a, described by the effective lagrangian 
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^ = + Ig'^'^d^ad^a, (3.113) 



levrG 2 

where we have added the Einstein-Hilbert action for the gravitational field as this also de 
scribes a low-energy degree of freedom. Note that with the potential defined in eqs.( |3.110| ), 



( p. 112 ) there is no cosmological constant in the presence of the vacuum value of the scalar 



field, p = py2/X or equivalently x = 0. 

The field equations derived from this lagrangian are the massless Klein-Gordon equation 
for the Goldstone field a: 

□"°V = -^d^./^g'"'d^a = 0. (3.114) 

and the Einstein-equations with a source term from the energy-momentum of the Goldstone 
field: 

R^^,. = -^^d^ad^a. (3.115) 

This system of equations admits exact planar wave solutions similar to those of the electro- 
magnetic waves discussed above. Using the same Ansatz (|3.54| ), (|3.76| ) for the metric we find 



running planar gravitational waves described by a metric component guu = —2K(u, y, z) 
in combination with a scalar wave (t{u) running in the same direction: 

□"°V = (-2 dud^ - 2d^Kd^ + dl + df) a = 0, (3.116) 

because both K and a are independent of v. The inhomogeneous Einstein equation now 
takes the form 



SttG 



Ruu = -AtransK = - ^^^^ {9^(7^, (3.117) 



C 



with the solution in complex transverse coordinates (C;C)- 



AttG 9 - - - 

K{u,y,z) = {d.crfa + f{u;0 + /(«;C), (3.118) 



c4 



where the (anti-)holomorphic terms Kq = f{u] () +c.c. denote the free gravitational waves 
that are allowed in the background. As the functional dependence of a{u) on u is quite 
general, we do not have to restrict ourselves to monochromatic waves, but can consider 
arbitrary wave packets. 

The net effect of a scalar wave on the motion of test particles is the same as that of an 
electro-magnetic wave. For the simple case of a monochromatic wave of wave number k: 

(r{u) = (To sin v^/cM, duCr{u) = kaoV2cos\/2ku, (3.119) 

the gravitational potential is of the form (|3.89|) with 
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^o'=^A;Vol (3.120) 

It follows that the solutions of the equations of motion for a point mass in a scalar wave is 
also given by expressions ( p.l02|) -( ^.104| ). 



In addition we also observe that scalar fields with zero frequency and constant gradient 
du<j{u) = 2^a/\fX represent a simple </)-wave 

0(m) = 00 e'""^ = e'°(-c*). (3.121) 

V A 

Such a wave can interpolate between two fixed values of the Goldstone field a on a finite 
interval —L < x — ct < 0: 

a{u) = (To, for u <0] 

/ N 2uau ^ ^ , r- 

a{u) = (To + forO < u < L/V2 



fiaLd—, foTU>L/\^. 



a[u) = o"! = ao 

The wave number a of the scalar field (f){u) is given by the slope of a in eg. ( p. 122 



as 



From the geodesic equation ( |3.96D with cji = we find, that in the scalar field with constant 
non-zero gradient — in the region < u < L/ a/2 — a point mass will execute a harmonic 
motion in response to the gravitational field generated by the uniform scalar potential 
difference. This response is the same as to the gravitational field of a uniform static 
electro-magnetic field. But here it is generated by a Goldstone field of constant growth, 
which really represents a constant change in the phase of a complex scalar field in the 
broken symmetry mode; as such the oscillations driving test masses through the associated 
gravitational field represent not only a mathematical, but also a physical analogue of the 
Josephson effect in superconductivity. 

From eqs.( |3.11"8D and ( p.78| ) it follows that these simple waves on the finite interval 
< u < L/\^ represent block-waves of constant curvature 

We have seen before that such block waves can be interpreted as dynamical gravitational 
domain boundaries; we observe that in this case they coincide with analogous domain 
boundaries of the scalar field. Stated differently, the transition between different scalar 
vacua also induces a transition between different gravitational vacua. 
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The simple treatment in terms of block waves can be regularized by considering more 
sophisticated field configurations interpolating between different vacua; for example the 
regular scalar field 

a{u) ^-{a, + a,) + ^^^^ arctan ( ^ , (3.125) 

2 TT \ V\ cTi - ao J 

interpolates between ctq and Ui; it also generates a smooth gravitational field: 



_1 _ 16tiG n'^a^ 1 

^ I 1 " 47r^M^ \ 

The geodesic equations describing the motion of test masses can be solved analytically 
in closed form for this space-time geometry. Observing that as before ii — ^ — const., 
the equations of motion in the transverse plane can be written conveniently in terms of u 
rather than proper time r as affine parameter. Then 



d?y UqU (Pz uJqZ 



where the parameters in the equation are given by the explicit expressions 



+ 

r 



, ^ sin ( Jl + Tj"^ arctan Tu] \ . 



(3.127) 



^0 = = w ^ 3.128 

We will also define the ratio of these quantities as a new parameter 

- ^. (3.129) 
We give the solution of the equation for y{u), that of z{u) being similar: 

y{u) — vT+T^t? (^Ho cos (^^J l-\-rf arctan Vu^ + 



(3.130) 



Here yo a-nd Vy represent the value of the ^/-coordinate and the velocity in the y-direction 
at = 0. Note that the velocity u — ±00 is given by 




(3.131) 
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Clearly the change in proper velocity depends on the position i/q: 

= 27ryocos (^^1 + ^') • (3-132) 

Prom this we see once more that a system of particles at fixed distances and at rest 
w.r.t. each other at some early time Ti — oo will be scattered by the wave and have 
all different velocities at a very late time r2 ^ oo: no Minkowski frame at early times 
remains a Minkowski frame at late times. This confirms our previous conclusion that 
the gravitational vacuum states, the Minkowski spaces before and after the wave, are 
inequivalent and that space-time can have a gravitational domain structure originating in 
the breaking of general co-ordinate invariance down to Poincare invariance, in addition to 
a domain structure originating in broken internal symmetry. 



Au,. 



dr 



dy 
dr 
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Chapter 4 
Black holes 



4.1 Horizons 

The example of gravitational shock waves/domain walls has shown that space-time in 
four dimensions is a truely dynamical arena for physics even in the absence of matter, 
although adding matter and other fields of force — for example in the form of scalar or 
electromagnetic fields — makes the possible structures even richer. In this chapter we 
turn to another class of very interesting structures in four-dimensional space-time: static 
or quasi-static solutions of the Einstein equations which represent gravitating extended 
bodies, whose static fields become strong enough to capture permanently anything that 
gets sufficiently close to the central core. As these quasi-static objects even capture light 
they are called black holes, but it has turned out that their properties are much more 
interesting and peculiar than this ominous but dull sounding name suggests. 

Black holes possess two characteristic features. The first, which has given them their 
name, is the existence of a horizon: a closed surface through which objects can fall without 
ever being able to return. The second is the presence of a singularity, a locus of points 
inside the horizon where the curvature becomes infinite. It has even been conjectured by 
some authors that these two properties are so closely linked that any space-time singularity 
should always be hidden behind a horizon. This conjecture bears the name of the Cosmic 
Censorship hypothesis. Although as an unqualified mathematical theorem the hypothesis 
is certainly not correct, it could be true for singularities that arise in realistic processes 
from the gravitational collapse of macroscopic bodies, like imploding massive stars. 

In this chapter we restrict our attention to the special black-hole solutions that might 
be called classical 'eternal' black holes, in the sense that at least outside the horizon 
they describe stationary gravitational fields, and do not require for the description of their 
behaviour the inclusion of their formation from collapsing matter. In fact they are solutions 
of the source-free Einstein or coupled Einstein-Maxwell equations and are characterised for 
an outside observer completely in terms of their mass, their angular momentum and their 
electro-magnetic charge. One might alternatively think of them as black holes that have 
been formed in the infinite past and settled in a stationary state in which all trace of their 
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previous history has been lost, e.g. by emitting gravitational radiation. 

Because of this restriction these lectures do not include a detailed account of the process 
of gravitational collapse as studied by astrophysicists, but they do allow one to isolate and 
study in detail the special aspects of the physics of black-holes. It may even be that these 
very special simple solutions are relevant to elementary particle physics near the Planck 
scale, but that leads one to consider the problems of quantum gravity, a subject we will 
not deal with at length in these lectures. The interested reader can find more information 
on many aspects of black-hole structure and formation in the literature ||35|, 



4.2 The Schwarzschild solution 

The simplest of all stationary black-hole solutions of the source-free Einstein equations is 
that for the static spherically symmetric space-time, asymptotically fiat at spatial infinity, 
first described in the literature by Schwarzschildf] We present a derivation of this 

solution making full use of the symmetries and their connection to Killing vectors explained 
in chapter |. 

In view of the spherical symmetry we introduce polar coordinates (r, 6', ip) in three- 
dimensional space and add a time coordinate t, measuring asymptotic Minkowski time at 
r ^ oo. The Schwarzschild solution is obtained by requiring that there exists at least one 
coordinate system parametrized like this in which the following two conditions hold: 
{{) the metric components are t- independent; 

{ii) the line-element is invariant under the three-dimensional rotation group SO (3) acting 
on three-vectors r in the standard (linear) way. 

The first requirement is of course equivalent to invariance of the metric under time-shifts 
t t + 6t. Thus formulated both conditions take the form of a symmetry requirement. 

In chapter (Q) we discussed how symmetries of the metric are related to Killing vectors. 
Rotations of three-vectors are generated by the differential operators L of orbital angular 
momentum: 



d n 9 ^ d n ■ 9 

Li = —smu?— — cot cos 09 — — , L2 = cos 09— — cot sm 09 — — , 
oO dip 06 dip 

(4.1) 

Similarly, time-shifts are generated by the operator 

L„ = |. (4.2) 

According to the discussion in sect.(|2.4D the coefficients of the symmetry operators = 
L^d^ are to be components of a set of Killing vectors La, A = 0, 1, 2, 3: 



^The solution was found independently at the same time in 
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D^Lau + D.LAf, = 0. (4.3) 
These can equivalently be characterised as a set of constants of geodesic flow as in eg. ( |2.45| ): 

Ja = L^P,. (4.4) 

The relation with the metric is, that the Poisson brackets of these constants of geodesic 
flow with the hamiltonian are to vanish: 

{Ja,H} = 0, H=]^g^^p^p,. (4.5) 

The general solution of these conditions for the hamiltonian H is t-independent, whilst it 
can depend on the angular coordinates only through the Casimir invariant 

L'=pI + -A-n- (4.6) 
sm d 

It follows that the metric must be of the form0 

ds^ = -h\r)de + g^{r)dr^ + e{r)dn^, (4.7) 

where dVl? = dO"^ + sin^ dip^ is the angular distance on the unit sphere, and the coefficients 
h{r), g{r) and k{r)) are functions of the radial coordinate r only. Indeed, the hamiltonian 
derived from this metric is 

1 2 1 2 

~ 2^^* ^ 2^^" ^ 2F' 

Note that the set of functions {h, g, k) can not be unique, as one can always perform a 
coordinate transformation r — > f(r) to another metric of the same general form ( ^.7| ) with 
coefficients {h,g,k). There are several standard options to remove this freedom; to begin 
with, we choose one in which the spherical symmetry is manifest and the spatial geometry 
is as close to that of flat space as possible. That is, we take r as the solution of 

k\r) = rV(^)- (4.9) 

Then the line element becomes 

ds^ = -h\r)de + g\r) (dr^ + r^d^f) = ~h^{r)dt^ + g\r)dr^ . (4.10) 

Clearly it is rotation invariant, and the three-dimensional space-like part of the line element 
is conformal to that of flat space in polar coordinates, with the conformal factor g{r) 



^In this chapter we employ natural units with c ~ 1. 
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depending only on the radius, and not on the orientation or the time. In the hterature 
the coordinates chosen here to parametrize the Schwarzschild metric are called isotropic 
coordinates. 

Up to the freedom of radial reparametrizations, the form of the line element (|4.10| ) is 



completely determined by the symmetry requirements. To fix the radial dependence of 
the remaining coefficients h{r) and g{r) we finally substitute the metric into the Einstein 
equations. For the Einstein tensor G^u = R^y — \/2g^yR we find that it is diagonal with 
components 

Gtt = ^ ^h[9 +-9) , 

\ r g ) 

Grr = (2h'g' + ^ + -ihgy], (4.11) 

gh \ 9 r J 

Gee = = -4 (h"9 + hg" + - {hg)' - ^] , 

sm U gh \ r g J 

In the absence of sources or a cosmological constant these expressions have to vanish. As 
boundary conditions at infinity, where space-time must become asymptotically fiat, we 
require that r —>■ oo implies /i — > 1 and g 1. Then we have three differential equations 
for two unknown functions, and therefore the equations must be degenerate. One can 
check, that this is the case if the following two relations are satisfied: 

hg = g + rg', [r'ihg)']' = 0. (4.12) 

Indeed, in this case the various components of Gf^y become proportional: 

9^ 1 1 

= ~/^^" " ^ ~ r2 sin^ e ^"^'^^^ 

The conditions ( [4.12| ) are met for any g{r) of the form 

gir) = a + - + 4, h{r) = (4.14) 

r g 



It turns out that eqs.( 4.l4 ) indeed provide solutions of the Einstein equations satisfying 
the correct boundary conditions. They are most easily obtained from the first Einstein 
equation (|4.11| ): for hg ^ it is possible to factor out all /i-dependence and obtain an 
equation for g only: 

99" + -99 - 1 9" = 0. (4.15) 
r 2 

There is a unique solution with the required normalization for r —>■ oo: 



9{r) = + ' (4-16) 
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where m is an undetermined constant of integration. Eq. ( [4.12[ ) then immediately provides 
us with the solution for h: 

h{r) = '^^—^. (4.17) 

^ ' 2r + m ^ ^ 

The complete solution for the static, spherically symmetric space-time in isotropic coordi- 
nates therefore is 

Clearly this solution matches Minkowski space at r — > cxd by construction, and h and g 
change monotonically (decreasing and increasing, respectively) as a function of r untill 
r = m/2. At that value h = 0, and the metric has a zero mode. Then the derivation of the 
solution we have presented breaks down, and its continuation to values r < m/2 requires a 
careful interpretation, discussed below. In any case, it turns out that this singularity of the 
metric is not physical, as invariants of the curvature remain finite and well-defined there. 
Indeed, in the following we construct coordinate systems in which the metric is perfectly 
regular at these points and can be continued to regions inside the surface r = m/2. In 
the isotropic coordinate system the apparent singularity of the metric for r = is another 
coordinate artifact. It disappears if one brings spatial infinity to a finite distance by an 
appropriate coordinate transformation. Details are given in the next section. 



4.3 Discussion 



Modulo radial reparametrizations r —>■ r(r) the line-element (|4.18|) is the most general 



static and spherically symmetric solution of the source-free Einstein equations matching flat 
Minkowski space-time at spatial infinity (r — >■ oo). Therefore in the non-relativistic weak- 
field range (w <C 1, r ^ m) it should reproduce Newton's law for the gravitational field of 
a spherically symmetric point mass of magnitude M, according to which the acceleration 
of a test particle at distance r is 

d^r GM , 
-TT = 4.19 



As the spherically symmetric solutions of the Einstein equations ( ^.18| ) are distinguished 



only by a parameter m, in the Newtonian limit this parameter must be related to the mass 
M, thereby representing the equivalent gravitational mass of the object described by the 
spherically symmetric metric ( 4.18|) . 



We first show that our solutions reproduce Newtons law in the limit r ^ m, w ^ 1 
(= c). The equation of motion for a test particle is the geodesic equation ( p.7|) . For the 
spatial coordinates this becomes 
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+ (terms ~ 0{v/c)) 

(4.20) 



Note that in the hmit w <C 1 one has dx'^ = dt ^ dr, whilst for r — > oo the factor h/g^ — > 1. 
Upon the identification m = GM we thus indeed reproduce Newton's equation ([4.19|) . 

Taking into account the dimensions of Newton's constant, as follows from ( |4.19| ) and 
is expressed in eg. ( |1-2| ), and our natural units in which c = 1, we find that m has the 
dimension of a length, consistent with the observation that m/r should be dimensionless. 
We have already noted that in isotropic coordinates the metric component h vanishes at 
half this distance: r = ro = m/2. Presently this only signals that the isotropic system 
cannot be extended beyond this range, but the surface r = tq does have a special physical 
significance: it defines the location of the horizon of a spherically symmetric black hole. 
This characteristic value ro = m/2 of the radial coordinate is known as the Schwarzschild 
radius. 



For the range m/2 < r < oo the Schwarzschild solution in the form ( |4.18|) is well-defined 



and it matches the asymptotic Minkowski metric for large values of r. We noted that on 
the Schwarzschild sphere r = m/2 the metric is singular, but as a solution of the Einstein 
equations it can be continued to smaller values of r and then it is again well-defined in the 
domain < r < m/2. However, it turns out that this domain of r- values does not in any 
sense represent the physical interior of the horizon, as one might expect naively. Rather, 
it easy to establish that the solution ( |4.18| ) is invariant under the transformation 



2 

m 



(4-21) 

which maps the domain (0,m/2) to (m/2,oo). Therefore the isotropic coordinate system 
actually represents a double cover of the region between the Schwarzschild radius and the 
asymptotic region r — > oo. Note that under this transformation h{r) —h{r), which can 
be compensated by an accompanying time-reversal t — > —t. From these observations it 
follows that the apparent singularity for r = disappears as it simply describes another 
asymptotic region like that for r — > oo in the original coordinate system. 

The existence of such a coordinate system presenting a double cover of the space outside 
the black hole might lead one to suspect, that a test particle thrown in radially towards the 
spherical Schwarzschild surface after a long enough time would return to an asymptotic 
Minkowski region. However, this does not happen; in the first place, as a simple calculation 
to be presented later shows, in terms of the time measured by an observer at asymptotic 
infinity it takes an infinite period of time for a test particle starting at any finite radial 
distance r > m/2 with any finite velocity even to reach the Schwarzschild surface; certainly 
such an observer will never see it start up and set off towards infinity again. Apart from 
this, leaving the Schwarzschild surface would also take an infinite amount of asymptotic 
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Minkowski-time. Clearly, these effects are due to the vanishing of h{r) on the horizon, 
which implies that an infinite amount of asymptotic Minkowski time t has to pass during 
any finite period of proper time r. But this is only the minor part of the argument; more 
important is the second reason, to wit that time-like geodesies do not flow from values 
r > m/2 to values r < m/2. Instead, as will become clear they flow into a new region of 
space-time described by complex values of r. 

Some calculational details are given later, but the above arguments at least indicate 
that an observer at spatial infinity is essentially disconnected from the domain of space- 
time beyond the Schwarzschild sphere. Therefore the Schwarzschild surface is called a 
horizon: from spatial infinity one cannot look beyond it. In contrast, things are radically 
different for an observer moving with the test particle towards the horizon. For such an 
observer in free fall, co-moving with the test particle on a radial geodesic and starting 
from an arbitrary point at a finite distance from the Schwarzschild surface, only a finite 
amount of proper time passes until the radius r = m/2 is reached, although for large radial 
distances the required time grows linearly with the distance. 

Moreover, having reached the horizon the particle (and the observer) can continue 
falling inwards without noticing anything particular, at least initially. The only remarkable 
effect they could discover once they have passed beyond the Schwarzschild surface is the 
impossibility to turn back to spatial infinity, no matter how powerful the engines they have 
at their disposal to accelerate and propel them. Therefore their passage to the inside of 
the Schwarzschild sphere represents an irreversible event: they are permanently trapped 
inside the Schwarzschild sphere. For this reason the Schwarzschild space-time is called a 
black hole: nothing, not even radiation, can get out once it has fallen through the horizon, 
at least not in the domain of classical relativity. 

To describe this motion of test particles (i.e. the flow of geodesies) through the Schwarzschild 
sphere, one must pass to a different coordinate system, one better suited to particles in free 
fall rather than to the description of the space-time from the point of view of an asymptotic 
Minkowskian observer at rest at spatial infinity. We make this passage in two steps. First 
we show that in other coordinate systems there is a different continuation of space-time 
beyond the Schwarzschild surface, one in which the interior geometry inside the horizon 
is distinct from the outside geometry and can not be mapped back isometrically to the 
exterior. After that we show by explicit calculation of the geodesies, that this alternative 
continuation can be taken to decribe the true physical situation for a test particle falling 
through the Schwarzschild sphere. 



4.4 The interior of the Schwarzschild sphere 

To go beyond the region of space-time covered by the isotropic coordinates, we perform 
a radial coordinate transformation which instead of making the spatial part of the metric 
conformal to flat space, as in (^4.9|) , ( 4.1(]| ), makes the angular part directly isomorphic to 
the two-sphere of radius f: 
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r = rg{r), k'^{r) = r^. (4.22) 

Then the metric takes the form 

ds^ = -¥{f)dt^ + g\f)df^ + f^dn\ (4.23) 
The exphcit form of the radial transformation ( [4.22| ) gives for the coefficients 

h\f) = ^ = 1 - ^. (4.24) 

This leads to the standard form of the Schwarzschild solution as presented in the original 
papers 

ds^ = - (^1 - ^) de + \^df^ + f^d^\ (4.25) 

f 



As to lowest order the transformation ( [4.22 ) does not include a rescaling, we have for large 
r and f: 

- = - + 0{\). (4.26) 

Therefore in the asymptotic region the Newtonian equation of motion for a testparticle in 
the field of a point mass M still holds in the same form after replacing r by f, and in the 
non-relativistic long-distance limit the two coordinate systems are indistinguishable. 

Note that the Schwarzschild radius r = m/2 is now located at f = 2m. Again the 
metric components are singular there, but now there is a zero mode in the (tt) -component 
and an infinity in the (rr)-component, in such a way that the determinant of the metric 
remains finite on the horizon. On the other hand, the determinant g does have a real zero 
for f = 0. This metric singularity was not present in the isotropic coordinate system, as 
it did not cover this region of space-time, but it turns out to represent a real physical 
singularity: the curvature becomes infinite there in a coordinate-independent way. 

A peculiarity of the metric ([4.23|) is that the components h"^ and g^ are not positive 
definite, but change sign at the Schwarzschild radius. This implies that the role of time 
and radial distance are interchanged for f < 2m: f is the time coordinate and t the 
radial coordinate. At the same time the metric components have become explicitly time 
dependent (they are functions of f) , and therefore the metric is no longer static inside the 
horizon. In return, the metric is now invariant under radial shifts t t + 6t. Hence it is 



the radial momentum pt which now commutes with the Hamiltonian (|4.8| ), the generator 
of proper-time translations, and is conserved. 

This strange result is the price one has to pay for continuing the Schwarzschild geometry 
to the interior of the Schwarzschild sphere. It is quite clear that this could never heve been 



achieved in the isotropic coordinate system ( 4.18 ), at least not for real values of r. Together 



with the appearance of the curvature singularity at r = this shows beyond doubt that 



65 



the region inside the Schwarzschild sphere r < 2m has no counter part in the standard 
isotropic coordinate system. 

To gain a better understanding of the difference between the continuation to r < m/2 in 
isotropic coordinates and f < 2m in Schwarzschild coordinates it is useful to perform yet 
another radial coordinate transformation to a dimensionless coordinate p defined by [^ | 



tanh^ = h{r), (4.27) 
2 

with the result that 

f = rg{r) = 2m cosh^ -. (4.28) 
Then the Schwarzschild metric becomes 

ds"^ = - tanh^ ^ dt^ + 4m2 cosh^ ^ (rfp^ + dQ'^) . (4.29) 

As coshp/2 > 1 this expression for the metric, like the isotropic coordinate system, is valid 
only outside the horizon f > 2m. Moreover, continuing to negative values of p we reobtain 
the double cover of the exterior of the Schwarzschild black hole: 

tanh(^) = -h{r) = h{^). (4.30) 

Now the interior of the Schwarzschild sphere can also be reparametrized in a similar way, 
taking account of the fact that h"^ has become negative: 

f = m(l + coscr) = 2m cos — , tan — = —h (f). (4-31) 

This parametrization clearly exists only for < f < 2m. In these coordinates the metric 
takes the form 

ds^ = tan^ I dt^ + 4m2 cos^ | [-da^ + dQ^) . (4.32) 

The time-like nature of a now is obvious. It is also clear, that this coordinate system 
presents an infinitely repeated periodic covering of the interior of the Schwarzschild sphere: 
the fundamental domain may be chosen to be < cr < vr, corresponding to 2m > f > 0; 
but continuation of a beyond this region leads to repeated covering of the same set of 
f-values. 

Comparing these two parametrizations of the exterior and the interior of the Schwarzschild 
sphere, it is now obvious that they are related by analytic continuation |29|: the interior 
of the sphere is not described by negative values of p, which correspond to r < m/2 in 
isotropic coordinates, but by imaginary values of p: 



a = ±ip, 



(4.33) 
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which is hke a Wick rotation changing the space-hke radial coordinate p to the time-hke 
coordinate a. In terms of the original isotropic coordinates it can be described as the 
analytic continuation of r to complex values 

Tfl 

r = -e'\ (4.34) 

Thus in the interior of the Schwarzschild surface the modulus of r is constant, but its 
phase changes, whilst outside its phase is constant and the modulus changes. Note that 
whereas the period of the fundamental domain of a in describing the interior region of the 
Schwarzschild sphere is vr, the values of the isotropic radial coordinate r are periodic in 
cr with period 2tt. Again it seems that the isotropic coordinates cover the Schwarzschild 
solution twice. We have more to say about that in the following. Finally we observe, that 



as indicated in eq.(|4.33|) there are two ways to perform the analytic continuation from p 



to o", i.e. from the exterior to the interior of the Schwarzschild space-time. Also this fact 
has non-trivial significance for the geometry of the Schwarzschild space-time. 



4.5 Geodesies 

The analysis presented so far makes it clear that the geometry of the Schwarzschild space- 
time is quite intricate. Indeed, at the Schwarzschild sphere r = m/2 (or f = 2m and 
o" = p = 0, respectively) several branches of the space-time meet: two branches of the 
exterior connected to an asymptotic Minkowski region, and two branches of the interior 
connecting the horizon with the space-time singularity. It is therefore of importance to 
understand the flow of geodesies in the neighborhood of the horizon, and see which branches 
test particles moving in the Schwarzschild space-time follow. We begin with time-like 
geodesies, representing the motion of test particles of unit mass. For such a particle one 
has in asymptotic Minkowski space: 

2H = g^^'p^p, = -1. (4.35) 



This relation was already argued from the definition of proper time in sect. (|2.1| ), eq. (|2.8| ) 
and below. As H is conserved, the relation holds everywhere on the particle's worldline. 

The solution of the geodesic equations for Schwarzschild space-time is much simplified 
by the large number of its symmetries. They provide us with sufficiently many additional 
constants of motion to allow a complete solution of the equations of motion of test particles 
in terms of their energy and angular momentum. The first of these constants of motion is 
the momentum component pt-, conserved because the metric is t-independent: 

,,^-h'(r)- = -. ^ ^^-j^y (4.36) 

At space-like infinity h'^ ^ 1; also, with r defined as proper time, dt = dr for a particle at 
rest there. Therefore a particle at rest at r — > oo has e = 1, equal to the (unit) rest mass of 
the test particle; moreover, as is intuitively obvious and argued more precisely below, if it 
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is not at rest it must have e > 1. Indeed, particles with e < 1 cannot reach spatial infinity 
and in this sense they are in bound states. The factor e gives the usual special relativistic 
(kinematical) time-dilation for a moving particle in asymptotic Minkowski space. 

In contrast, for finite radial distance vq < r < oo the factor /i^(r) describes the grav- 
itational redshift: the time dilation due only to the presence of a gravitational field. It 
represents a universal effect, not influenced by the state of motion of the particle. 

The conservation of angular momentum gives us three conservation laws, which are 
however not independent: 

Jl = — sin (f pq ~ cot 9 cos (f p^, 

J2 = cosippe — cot9 simpp^, (4.37) 

J3 = Pv" 

with Jl + J2 + Jl = L"^ = (?i a constant. Here the separate momentum components are 

pe = k\r)^, p^ = k\r) sin' 9^. (4.38) 

If one orients the coordinate system such that the angular momentum is in the x^-direction, 
then 9 = Ti /2 = constant, and therefore pe = Q and Pip = J^ = L Then the angular 
coordinates as a function of proper time are given by 

It remains to solve for the radial coordinate^] r. For this we use the conservation of the 
hamiltonian, setting 2H = — 1: 



Note, that consistency of this equation in the region outside the horizon requires 

^ > II + 1. (4.41) 

It shows, why in the region r — 00 the time dilation factor has to be greater or equal 
to unity: e > 1. At the same time, for non-zero values of angular momentum C. the 
radial distance r can in general not become arbitrary low: like in the ordinary Kepler 
problem there is a centrifugal barrier to be overcome. For details we refer to the literature 
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Combining eq.( [4.40| ) with that for the angle v?, by eliminating r we directly obtain the 
orbital equation 



^Note that each of these equations is still true for arbitrary radial parametrization. 
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It is clear that an orbit can have a point of closest approach, or furthest distance, only 
if dr/dip = dr/dr = 0. The only physical solutions are points where relation ( [4.41|) 
becomes an equality. When the equality holds identically, for all points of the orbit, the 
orbit is circular. One can show that the smallest radius allowed for a stable circular orbit 
corresponds in Schwarzschild coordinates to f = 2miQ = 6m. 

To the set of time-like trajectories we add the light-like geodesies, representing the orbits 
of massless particles. In that case the proper-time interval vanishes identically on any 
geodesic, that is to say H = 0. We parametrize the trajectories with an affine parameter 
A and obtain 

P. = -/.^^--7, (4.43) 
allowing one to replace the affine parameter by t through 

dX = — dt. (4.44) 
7 

One can again choose the plane of motion to be ^ = 7r/2, and the conservation of angular 
momentum holds in the form 

7r^ dip o dip .rN 
where is a constant. The radial equation now becomes 

or using the choice 6 = n/2 above: 

Like in the massive case, if 7^ one can eliminate the affine parameter A also in favour 
of p) instead of t. 

An interesting case, as concerns the exploration of the horizon and the connection between 
the interior and the exterior of the Schwarzschild sphere, is that of radial motion: £ = or 
u; = 0, respectively. Let us first consider the time-like geodesic for a particle starting from 
rest near r = 00, which implies e = 1 and dr/dt < 0. Then 

dt 1 dr 1 
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It is easy to solve these equations in the Schwarzschild coordinates; indeed, with {l — h"^) = 
2m /f one finds a solution valid in the whole range < f < oo, both outside and inside the 
horizon: 



2m 



4m 



{ro - r) 



(4.49) 



where Tq is the proper time at which the particle reaches the singularity at f = 0. Taking 
this time as fixed, r denotes the proper time, prior to tq, at which the particle is at radial 
distance f. Observe, that this time is finite for any finite radial distance f; only when 
starting from infinity it takes an infinite proper time to reach the singularity. Nothing 
special happens at the horizon f = 2m: it is crossed at proper time 



4m 



(4.50) 



a proper period 4m/3 before the singularity is reached (of course, in terms of the time- 
coordinate f inside the horizon, it takes a finite period 2m). 

Solving for i as a function of proper time, we have to distinguish two cases: f > 2m 
and f < 2m. Namely, introducing the new variable 



y 



we obtain the differential equation 



4m 



1^0 



(4.51) 



1 dt 
4m dr 



y 



(4.52) 



the solution of which depends on whether the r.h.s. is positive or negative. We find for the 
two cases the solutions 



4m 



Tq-t 
4m 



+ 



4m 



{ro - r) 



+ 



arccoth 
arctanh 



4m 
3 

.4m 



{ro - r) 
{ro - r) 



r > 2m; 
r < 2m. 



(4.53) 



In both cases the left-hand side becomes infinite at the horizon, i.e. at proper time th- 
This confirms that an external observer at infinity will never actually see the particle reach 
the horizon. Similarly, the distance measured from the horizon to the singularity at r = 
is infinite in terms of the coordinate t. 

Now consider what happens in isotropic coordinates. In terms of a new variable 



(4.54) 
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where we choose the positive branch of the square root, integration of the equation for r(r) 
is elementary and gives 

l.^r„-r)^l(. + l^f. (4.55) 

Clearly this is symmetric under interchange x — corresponding to inversion w.r.t. the 
horizon r = m/2. Indeed, both at r — cxd and r = we see that r — > — oo. This confirms 
our earlier remark that a particle released from space-like infinity does not return there. 
Instead, it will cross the horizon x = 1 at r = r//, as in the previous calculation. Then 
continuing r to complex values (|4.34| ), we find 



x = e5^ !^=icos3-. (4.56) 

4m 3 2 ^ ' 

The singularity a = vr is then reached for r = Tq, as expected. In terms of the radial 



coordinate p, related by the pseudo-Wick rotation (|4.33|) the equation for the radial distance 



in terms of proper time outside the horizon becomes 

I^ = ^cosh^P. (4.57) 
4m 3 2 ^ ' 

The analytic continuation of coordinates in the complex plane on crossing the horizon is 
manifest . 

Next we consider incoming light-like radial geodesies. The equation is 

dr h , ^ 

- = — . (4.58) 
dt g 

Let to refer to the time at which the massless particle (e.g., a photon) passes a fixed point 
on the radius, e.g. fo in Schwarzschild coordinates. Then the solution of the geodesic 
equation for t is 



t-to ro-r 

= — m 

2m 2m 



2m 



ro — 2m 



(4.59) 



With the absolute value as argument of the logarithm this equation is valid both outside 
and inside the horizon. The time for the photon to reach the horizon from any point outside 
diverges, as does the distance it has to move from horizon to the curvature singularity, as 
measured by the coordinate t. 



4.6 Extended Schwarzschild geometry 

As we have already noticed that for a massive particle falling across the horizon nothing 
very special happens, it is unsatisfactory that we do not have a non-singular description of 
the same process for massless particles. To improve the description of light-like geodesies at 
the end of the last section, it is not sufficient to perform a radial coordinate transformation: 
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one has to ehminate the t-coordinate and replace it by one which does not become infinite 
on the horizon. Moreover, to describe hght-hke geodesies it would be convenient to have 
a metric where the light-cone structure is manifest by having the radial part of the metric 
conformal to radial Minkowski space: a coordinate system in which (r, t) are replaced by 
coordinates {u, v) such that the radial part of the metric becomes 



(4.60) 



Such a coordinate system was first proposed by Kruskal and Szekeres [ED|; defining the 
new coordinates such as to be dimensionless, the line element is written in the form 



ds^ 
Am? 



f{u,v) (^-dv'^ + du^^ + g{u,v)dVt^ 



(4.61) 



with the standard choice for (m, v) defined in terms of the Schwarzschild coordinates (r, t) 
by: 



2 2 
U — V 



{ — 1 ) 62"- , 

V2m ' 



u 

V 



tanh — , if If I < \u\ 
Am 



coth - — , if If I > \u\. 
Am 



(4.62) 



In these coordinates the horizon f = 2m, t = ±oo is located at m = ±f , whilst the curvature 
singularity f = corresponds to the two branches of the hyperbola v? — v"^ = —1, one in 
the past and one in the future of an observer at space-like infinity. More generally, from 
the first result one concludes that the surfaces f = constant are mapped to hyperbola's 

— v"^ = constant, where the last constant is positive outside and negative inside the 
horizon. These hyperbola's are worldlines of particles in circular orbit, i.e. being subject 
to a constant acceleration. 

Similarly hypersurfaces t = constant correspond to u = constant x v, represented by 
straight lines in a {u, f )-diagram. Note also that the singularity f = consists of a set of 
two space-like surfaces in {u, f )-coordinates. 

More explicitly, we can write {u, v) in terms of (r, t) as 



f j_ , t 

u = ±\ le4m cosh- — . 

V 2m 4m 



(4.63) 



f j_ , t 

V = ±W le4m smh , 

V 2m 4m 



outside the horizon (r > 2m); and 
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/ r j_ t 

u = ±\/l 64"! sinn - — , 

V 2m 4m 



(4.64) 



f j_ t 

±W 1 64"! cosh- — , 

V 2m 4m 



within the horizon (f < 2m). We note once again, that these coordinates produce two 
solutions for each region of the space-time geometry. In the above conventions the dimen- 
sionless coefficient functions are given imphcitly by 



f[u,v) = —e g{u,v) = (4.65) 



We can also compare with the coordinate systems ( [4.29| ) and ( [4.32[ ), describing explicitly 
the double cover of the Schwarzschild space-time: 

u = cosh2 £ j^j^ P _^ V = ^ sinh - sinh — , (4.66) 

2 4m 2 4m ^ ^ 

for the exterior region \v\ < \u\; and 

u = 2 sm — smh — , f = 2 sm — cosh , (4.67) 

2 4m 2 4m ^ ^ 

for the interior region \v\ > \u\. From these expressions we infer, that the exterior regions 
p > and p < 0, corresponding to r > m/2 and r < m/2 in isotropic coordinates, and 
the two interior regions a > and a < as well, coincide with the regions [u, v) > 
and {u, f ) < in Kruskal-Szekeres coordinates, respectively. Hence the double cover of 
Schwarzschild space-time given by the two solutions for the Kruskal-Szekeres coordinates 
above correspond exactly to the double cover we found in isotropic and (p, cr)-coordinates. 

The Kruskal-Szekeres coordinates were introduced, and are especially useful, because 
they represent the light-like geodesies in a very simple way. In particular, the radial light- 
like geodesic are given by 

du^ = dv^ ^ = ±1. (4.68) 

dv 

These correspond to straight lines in the {u, f )-diagram parallel to the diagonals. The 
special light rays u = ±v are the asymptotes of the singularity and describe the horizon 
itself: they can be interpreted as photons permanently hovering on the horizon. Any other 
light-like radial geodesic eventually hits the singularity, either in the past or in the future, 
at finite [u, v), although this is in the region before t = —00 or after t = cxd in terms of the 
time of an observer at spatial infinity. 

As the light cones have a particularly simple representation in the {u, f )-diagram, for all 
regions inside and outside the horizon, it follows that time-like and space-like geodesies are 
distinguished by having their derivatives \du/dv\ > 1 for the time-like case, and \du/dv\ < 1 
for the space-like ones. One can check by explicit computation, that a particle starting from 
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rest at any distance r and falling into the black hole has a time-like world line everywhere, 
also inside the horizon: du/dv can change on the world line of a particle, but never from 
time-like to space-like or vice versa. Therefore all time-like geodesies must cross the horizon 
(light-like) and the singularity (space-like) sooner or later. If they don't reach or start from 
spatial infinity, they must come from the past singularity across the past horizon, and then 
cross the future horizon to travel towards the final singularity, all within a finite amount 
of proper time. 

The Kruskal-Szekeres coordinate system also teaches us more about the topology of the 
Schwarzschild geometry: how the various regions of the space-time are connected. First of 
all it shows that there is a past and a future singularity inside a past and a future horizon 
(|m| < |f I), depending on v being positive or negative. In this respect it is helpful to observe, 
that according to eq. (|4.64|) for f < 2m the sign of v is fixed for all t. In the literature these 



regions are sometimes refered to as the black hole, the interior of the horizon containting 
the future singularity, and the white hole containing the past singularity. In a standard 
convention, these regions of space-time are classified as regions II (future Schwarzschild 
sphere) and IV (past Schwarzschild sphere), respectively. 

There are also two exterior regions connected to asymptotic Minkowski space. They 
are distinguished similarly by m > and m < 0, as for f > 2m the sign of u is fixed for 
all t: one cannot reach negative u starting from positive u or vice versa by any time-like 
or light-like geodesic. In the same standard convention these regions are refered to as 
region I and region III, respectively. Being connected by space-like geodesies only, the two 
sheets of the exterior of the Schwarzschild space-time are not in causal contact and cannot 
physically communicate. However, the surfaces t = constant do extend from m = cxd to 
u = —oo] indeed, the regions are space-like connected at the locus u = v = 0, where 
the two branches of the horizon intersect. To understand the topology of Schwarzschild 
space-time, we must analyze this connection in some more detail; this can be done using 



the technique of embedding surfaces, as explained e.g. in ref.|TT| . 

Consider a surface t = constant and fix a plane 6 = 7r/2. We wish to understand how 
the regions u > and u < are connected for this plane. The Schwarzschild line element 
restricted to this plane is 

ds' = + f^V. (4.69) 

f 

Now we embed this two-dimensional surface in a three-dimensional fiat, euclidean space: 

ds^ = dx^ + dy^ + dz'^, (4.70) 
by taking x = fcosip, y = fsimp, and defining the surface z = F{f) such that 

ds^ = (l + F'\f)) df^ + f^d^\ (4.71) 
Comparison of eqs. (|4.71|) and ( [4. 691) then implies 
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16m"' y 4m J 2m 

This defines a surface of revolution, obtained by taking the parabola z"^ = 8mx — 16m^ in 
the {x, z)-plane and rotating it around the z-axis. It is a single surface, with looks like a 
throat, of which the upper half z > corresponds to 





or M > 0, (4.73) 
whilst the lower half z < describes the branch 

or M < 0, (4.74) 

The two regions are connected along the circle f = 2m at 2; = 0. Of course, the spherical 
symmetry implies that the same topology would be observed for any choice of angle 6; 
any plane with -u > is connected to one with m < in this way. Although no time-like 
travel is possible classically between the regions on each side of the throat, the possibility 
of particles tunneling via quantum processes between the two regions poses a curious and 
interesting question to theories of quantum gravity. 



4.7 Charged black holes 

The spherically symmetric solution of the source-free Einstein equations can be extended 
to a solution of the coupled Einstein-Maxwell equations including a static spherically sym- 
metric electric field (in the rest frame of the black hole)[^. As the only solution of this 
type in flat Minkowski space is the Coulomb field, we expect a solution of the coupled 
equations to approach the Coulomb potential at space-like infinity: 

A = A.dx" = ^rft, (4.75) 

4:Tcrp[r) 

with p{r) — * 1 for r 00. This is a solution of the Maxwell equation 

D^F^"" = 0, (4.76) 

in a spherically symmetric static space-time geometry. The Maxwell field strength appear- 
ing in this equations is 



F.^ dx^ Adx" = ^ ^'T ; dt A dr. (4.77) 



2 47r r^p^ 

In isotropic coordinates ( 4.1CI|) the energy- momentum tensor of this electro- magnetic field 
reads 
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T,AF] = ( F,,F/ - i g,.FA = ^^±# diag -1, r^ W ^"l . (4.78) 



This has the same general structure as the Einstein tensor ( [4.13| ): 



G,. = ^ Uh'g' + ^ + l (hg)'] diag -1, r^ sin' 9 ] , (4.79) 



which was derived under the conditions ( [4.12| ), with the general solution ( [4.14| ). Now 
requiring the coupled Einstein-Maxwell equations 

G^, = -8nGT^,[F], (4.80) 
this reduces to a single condition on the radial functions {h,g,p), with the solution 

2 

771 V 

P{r) = g{r) = 1 + — + — , 

(4.81) 

h(r) = = - 

g 4r^ + 4mr + u'^ 

Here the constant z/^ is defined by 

2 2 2 2 y '-^ /A Qn\ 

u = 7n — e , e = . (4.82) 

4:71 

It may be checked that this solution satisfies the Maxwell equations in this gravitational 
fields as well. As g{r) can also be written in the form 

we find that the Schwarzschild solution is reobtained in the limit e' 0. For z/^ > the 
zero of the metric coefficient h'^ is now shifted to 

r = ^. (4.84) 

As in the Schwarzschild case this is a point where two asymptotically fiat regions of space- 
time meet: the metric is again invariant under the transformation 

r - ^, t^-t, (4.85) 

under which h{r)dt — > h{r)dt. Moreover, also like in the Schwarzschild geometry, the 
sphere r = z//2 actually defines a horizon and time-like geodesies are continued into the 
interior of this horizon, which is not covered by the (real) isotropic coordinate system. 
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A Schwarzschild-like coordinate system which allows continuation into the interior of 



the horizon exists. It was found by Reissner and Nordstrom and is related to the 
isotropic coordinates by a radial coordinate transformation 

r = rg(r) = r + m + —. (4.86) 
Then the metric of the charged spherically symmetric black hole takes the form 

ds'^ = -h^dt^ + fdf'^ + f^dn^ (4.87) 
where the coefficients {h, g) are again inversely related: 

The horizon is now seen to be located at r = m + z/; in fact, there is a second horizon at 
f = m — u, provided m > e. As in general h"^ and g'^ are not positive definite, at each of 
these horizons the role of radial and time coordinate is interchanged. Only for u = 0, i.e. 
for m = e, the situation is different, as the two horizons coincide and no interchange of 
time- and space-like coordinates occurs. This special case of Reissner-Nordstrom solution 
is refered to in the literature as an extremal charged black hole. 

The result of the second horizon for non-extremal charged black holes is, that the real 
curvature singularity that is found to exist at r = is a time-like singularity, for m > e. 
This is expected from the Coulomb solution in fiat Minkowski space, which contains a 
time-like singularity of the electric field at r = 0. Therefore such a black hole looks more 
like a point particle than the Schwarzschild type of solution of the Einstein equations. 

The interior region between the horizons can be obtained from the isotropic coordinates 
by complex analytic continuation of r from the outer horizon by 

r = ^e'". (4.89) 

The domain of the argument a is chosen as [0, tt]. Indeed, for these values of r the radial 
Reissner-Nordstrom coordinate is 

f = m + ucosa, (4.90) 

and takes all values between the two horizons f = m±z/ precisely once. In this parametriza- 
tion the isotropic Reissner-Nordstrom metric becomes 



2 • 2 

ds"^ = — ^ ^ ^ dt"^ + (m + ucosaf i-da^ + dn'^) . (4.91) 
im + z/coscr) ^ ^ 



Again we observe the role of a as the time-like variable, with t space-like between the 
horizons. For o" = vr the isotropic coordinate r becomes negative: 

a = TT ^ r = (4.92) 
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Thus we discover that the Reissner-Nordstrom space-time can be extended consistently to 
negative values of the radial coordinate. In particular the curvature singularity at f = is 
seen to be mapped to 

^ ' (4.93) 



2 

This shows first of all, that there are two (time-like) branches of the singularity, and 
secondly that r can become negative for non-extremal black holes with m > e. The two 
branches are found also upon analytic continuation of a in ( [4.91|) to imaginary values: 



a ±ip. (4.94) 

This leads to a parametrization of the Reissner-Nordstrom geometry for the exterior regions 
of the form 

2 . 1 2 

ds^ = ^ 2 + (m± i/coshp)' (dp^ + dQ^) , (4.95) 

(m ± cosh p) ^ ^ 

where the plus sign holds for the exterior region connected to asymptotic fiat Minkowski 
space, and the minus sign to the region connected to the curvature singularity. Because 



of the double- valuedness of the analytic continuation ( [4. 941) there are two copies of each of 
these exterior regions involved. 



4.8 Spinning black holes 

The Schwarzschild and Reissner-Nordstrom black holes are static and spherically symmetric 
w.r.t. an asymptotic Minkowski frame at spatial infinity. There also exist extensions of 
these solutions of Einstein and coupled Einstein-Maxwell equations which are neither static 



nor spherically symmetric, but represent rotating black holes with or without charge [42 



43| . These solutions have an axial symmetry around the axis of rotation, and are stationary: 
in the exterior region the metric is time- independent, and the rate of rotation and all other 
observable properties of the black hole are constant in time. 

The standard choice of coordinates for these metrics, with the axis of rotation being 
taken as the 2;-axis and reducing to those of Schwarzschild, or Reissner and Nordstrom for 
the spherically symmetric case of no rotation, is that of Boyer and Lindquist |Q . Denoting 



these coordinates by (t, r, 6, if) to emphasize their relation with the spherical coordinate 
systems for the non-rotating metrics, the line element for reads 



ds^ = — — (dt- asm^9d(p) + ^ (df^ + A^dO' 

r'' + a'' cos^ 9 ^ ' ^ 

(f2 + a^f sin^ 9 ( ^ a n ^ 

+ -2 I 2 ^ dip - — — ^dt 

+ cos"' 9 \ r-^ + a-^ 



(4.96) 
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Here a is a constant parametrizing the deviation of this hne element from the usual diagonal 
form of the Schwarzschild-Reissner-Nordstrom metric; its physical interpretation is that of 
the total angular momentum per unit of mass: 

J = ma. (4.97) 

As one can always choose coordinates such that the angular momentum points along the 
positive 2;- axis, the angular momentum per unit of mass may be taken positive: a > 0. For 
convenience of discussion this will be assumed in the following. 
The quantity in the line element ( [4.96| ) is short hand for 

= {f-mf + + - ml (4.98) 

Note that, like the metric coefficients h"^ and we have encountered earlier, this quantity 
is not positive definite everywhere. Also, the zero's of this function determine the location 
of the horizons, of which there are two (an outer and an inner one) for > + e^: 



r±H = m ± \Jm^ — — e^. (4.99) 
The gravitational field represented by the metric ( [4.96|) is accompanied by an electro- 



magnetic field, the vector potential of which can be chosen as the one-form 

A = (dt - a sin^ Odip) . (4.100) 

For convenience of notation we have here introduced the function p defined as 

f = + (4.101) 



These fields are solutions of the coupled Einstein-Maxwell equations ( [4.76|) and (|4.80| ), with 
the Maxwell field strength given by the two form 

^f^ — cos^ 6^ dr A (^dt — a sin^ 



47rp4 



(4.102) 

qa r cos 9 sin 6 



+ I dO A (adt- (H + a^)dip) 



2tt p 

and with the identification of the quantity as in eq.( [4.82D : 



Air ' 



To determine the properties of the space-time described by the Kerr-Newman line element 
we need to analyse the geodesic fiow. Because the spherical symmetry of the Schwarzschild 
and Reissner- Nordstrom solutions is absent, and only axial symmetry remains, there is one 
less independent constant of motion from angular momentum: only the component along 
the axis of rotation (the 2;-axis, say) is preserved. Fortunately this is compensated by the 
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appearance of a new constant of motion, quadratic in momenta, which replaces the Casimir 
invariant of total angular momentum [^, As a result one can still completely solve 
the geodesic equations in terms of first integrals of motion. We now give some details. 
First we write down the geodesic hamiltonian, effectively representing the inverse metric 



2 

(4.103) 



Clearly, to be well-defined we must assume 6 ^ (0, tt), as the z-axis represents a coordinate 
singularity in the Boyer-Lindquist system where the angle ip is not well-defined. The above 
hamiltonian is a constant of geodesic flow; recall, that for time-like geodesies it takes the 
value 2H = —1, whilst for light-like ones H = 0. 

Next there are two Killing vectors, respresenting invariance of the geometry under 
time-shifts and rotations around the z-axis: 



-Pt - -9ttYx-3t^dX 



dt dp 

(4.104) 



dt dip 
I = P, = + 9,,^. 

Of course, for time-like geodesies the affine parameter A maybe taken to be proper time 
T. These equations can be viewed as establishing a linear relation between between the 
constants of motion {e,£) and the velocities {dt/dX,dp/dX): 

(4.105) 

with the dot denoting a derivative w.r.t. the affine parameter A. The determinant of this 
linear transformation is 

detG = -gug^^ + g^ = sin^^, (4.106) 

and changes sign whenever A^ does. Therefore the zero's of A^ define the locus of points 
where one of the eigenvalues of G changes sign, from a time-like to a space-like variable 
or vice-versa; at the same time, the coefficient of in the hamiltonian ( [4.103| ) changes 




sign, confirming that these values represent the horizons of the Kerr-Newman space-time. 
It follows from this analysis, that outside the outer horizon, in the region connected to 
space-like infinity, det G > 0. 

Now if from eqs. (|4.104|) the angular proper velocity dip/dX is eliminated, one finds the 
relation 

^„ detG dt 
e - Vl£ = — , (4.107) 

g^^p dX 
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where the quantity with the dimensions of an angular velocity is defined by 

f2 + - 
^ (f2 + a2)2 - sin^^ 

(4.108) 

1 

1 + A2p2 (f2 + ^2)-! (^2 + ^2 _ A2)-l " 

Observe, that it is proportional to a with a constant of proportionality which for large but 
finite f is positive. Hence Vl represents an angular velocity with the same orientation as 
the angular momentum J. With a canonical choice of coordinate system this quantity is 
positive, and therefore > as well. 

It also follows from eg. ( [4.107|) that for time-like geodesies outside the horizon, with 
A = r and therefore dt/dX > 0, an inequality holds of the form 



n 



9w 



r2 + a? 



e-fi£>0. (4.109) 

For large f near infinity one finds in fact the stronger inequality e — Vtl > 1. 

Returning to the solution of the geodesic equations, we take into account the existence of 
a constant of motion K such that 2K = K^'^p^p^, with K^" the contravariant components 
of a Killing tensor: 



Df^^K^x) = 0. (4.110) 

The explicit expression for K is 



1/a22 2/>2 2 2\ ^'^ sin^ 6 f p^p 



2p2 

cos^ 9 



2-p 



7.2 



apt 



(4.111) 



+ 



2p2A2 



((r^ + a^)pt + ap^ 



That K defines a constant of geodesic fiow is most easily established by checking that its 
Poisson bracket with the hamiltonian H (|4.103|) vanishes: 



{K,H} = 0. (4.112) 

Using these constants of motion, it is straightforward to establish the solutions for time- like 
geodesies in the form 
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dt 
dr 



A2 



+ e- 



(f2 + a2)2 



A^a^ sin^ 6 



A2 



_2dip 



' dr^ 
dr 



A^ - sin^ 6* 
A2 sin^ ^ 



— 1 — 



A2 



-A^ (2fs: + f2) + (fa - £(f2 + 



(4.113) 



2fs: - cos^ ^ - sin^ ^ 



a sin^ 



Note that for a particle starting from r = cxd one has £ > 1. These solutions are well-defined 
everywhere except on the horizons and at points where = 0. At the latter locus the 
curvature invariant Rp^ppaR^^'^" becomes infinite and there is a real physical singularity. 
At the horizons both t and (p change infinitely fast, signalling infinite growth of these 
coordinates w.r.t. an observer in asymptotic flat Minkowski space. 

As with the infinite amount of coordinate time t to reach the horizon in the Schwarzschild 
and Reissner-Nordstrom solutions, this merely signifies a coordinate singularity. However, 
here it is not only t but also (p which becomes infinite. To gain a better understanding of 
the reason for this rotation of geodesies, first notice that according to the second equation 
( [4.113| ) even for i = the proper angular velocity dip/dr does not vanish in general, be- 
cause e 0. Hence there are no purely radial geodesies: any incoming test particle picks 
up a rotation. Now computing the angular velocity as measured by an observer at rest at 
spatial infinity from the first two equations ([4.113|) one finds: 



dip 
'dt 



£9w + ^9t<p 

It follows immediately, that for a geodesic without orbital angular momentum: 
observed angular velocity equals the special reference value (|4.108| ): 



dip 
~dt 



At the horizon r = r^n, where A^ = 0, this becomes 



dip 
'dt 



£ = 

f = f+H 



H 



-H 



+ a 



2 ■ 



(4.114) 
= 0, the 

(4.115) 
(4.116) 



Comparison with the expression for Q, eq.( [4.108[) , shows that outside the horizon the 
inequahty 



n < fin, 



(4.117) 
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holds, with equahty only on the horizon. As these quantities are finite, the time t can 
increase near the horizon at an infinite rate only if ip increases at a proportional rate. For 
this reason both t and ^9 have to grow without bound as f approaches f+//. 

At space-like infinity the Killing vectors ( |4.104| ) define time translations and plane 
rotations. However, their nature changes if one approaches the horizon. Consider the time 
translation generated by pt = The corresponding Killing vector has contravariant and 
covariant components, given by 

= -S^ ^ 6, = -9t,. (4.118) 

Here 6^ is the usual Kronecker delta symbol. Similarly the components of the Killing 
vector for plane rotations, generated by p^, are 

^!; = si; ^ u = 9^,- (4-119) 

It follows, that these Killing vectors have norms given by 

{it? = 9tt = , 4.120) 

P 



and 



(U = 9^^ = sin' 6^- ^-^ , (4.121) 

whilst their inner product is 

^f^^ = -9t^ = asm'e^-^^^-^^. (4.122) 

From these equations several observations follow. First, as the norms and inner products of 
Killing vectors are coordinate independent quantities, we can write the expressions for the 
reference angular velocity fl, eq.( [4.108| ), and for the observed angular velocity, eq.( |4.114] ), 
in a coordinate independent form: 

a ' dt- ile-ifi,f ^'-'"'^ 

The last equality is equivalent to 

(.-Jifl^ = n.-.| = n(.-n^^). (4.124, 

Second, with a > and excluding the 2;-axis 9 = (0, vr) — where the Boyer-Lindquist 
coordinates are not well-defined and Q vanishes — Q and the invariant inner product C,t ■ 
C,ip are positive outside the horizon. Also the positivity of det G in the exterior region, 
eq. ([4.106|) , can be expressed in the invariant form 
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^'>^ ^ ^<^^ (4-125) 



with the equahty holding only on the horizon. Then for £ > we find from eg. ( [4. 12^ ) 



that dip/dt > whilst for i < the angular velocity is less than this reference value: 
d(f/dt < Q. As long as is time-like this is all that can be established; in fact, with i < 
and < 0, the absolute value of the angular velocity dip/dt can become arbitrarily large. 
However, it is easily established that the norm of C,t changes sign on the surface 



f{e) = fo{e) = m + Vm2-e2-a2 cos^^, (4.126) 

which, with the exception of the points where the surface cuts the z-axis: 9 = (0,7r), 
lies outside the horizon. The region between this surface and the horizon is called the 
ergosphere, and C,t is space-like there: > 0. It then follows directly from eqs.( f4.12"^ and 
( [4.1251 ) that in this region the angular velocity is always positive: 

^ > 0, for > 0. (4.127) 

Thus we have established, that inside the ergosphere any test particle is dragged along 
with the rotation of the black hole and no free-falling observer can be non-rotating there. 
Related to this, both and C,^ are space-like Killing vectors there. Nevertheless, outside the 
horizon there still exists time-like combinations of these vectors even inside the ergosphere: 

X = 6 - ^e^, (4.128) 



which has non-positive norm: < as a result of the inequality ( [4.125| ), with equality 
holding only on the horizon f = fn- 



4.9 The Kerr singularity 

We continue the analysis of black-hole geometry with a brief discussion of the nature of 
the singularity of spinning black holes. To begin, we have already noticed that, like the 
spherical charged black hole, the spinning Kerr-Newman black holes have two horizons, an 
innner and an outer one. As the singularity = is located inside the inner horizon, we 
expect the singularity to have a time-like character. 

In order to elucidate the topological structure of the space-time near the singularity, 
we introduce new coordinates in two steps. First we define the Kerr coordinates (V^, (p) as 
the solution of the differential equations 

dV = dt + to^ dr, d(p = d^ + -^df. (4.129) 

A A 

With these definitions the Kerr-Newman line element becomes 
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ds^ = — — (dV - asi-n^ edii) H — \if + a^)di2)-adV 



(4.130) 



+ 2 (dt/ - a sin^ Q d(p) df + f- dd'^ . 



Now we turn from quasi-spherical coordinates to quasi-Cartesian ones, known as Kerr- 
Schild coordinates, by defining 

X — [f cos — a sin (p) sin ^, 

y — (f sin -I- a cos sin ^ (4.131) 

z = rcos9, 

and we reparametrize the time coordinate t ^ t using 

£ = dV - df. (4.132) 
We then find that at fixed {9, if) the transformation 

X x' — X — 2fcos<^, 

y —>■ y' = y — 2fsin(^, (4.133) 
z ^ z' — — 

is equivalent to f — >■ f' = — f. This shows, that like the Reissner-Nordstrom case, the 
negative values of f are to be taken as part of the physical domain of values. 
Next one introduces the coefficients (i = 1, 2, 3) by the relation 

k-dv^ -ri^dx + ydy) + a{xdy-ydx) _ z_dz _ 

-\- r 



and scalar function u{r): 



, , -A^-Fr^-Fa^ 2mr - /.^o.x 
u{r) = ^ = (4.135) 

J.Z _|_ 



With these definitions the metric can be written in the hybrid form 

ds^ = dxl + u{r) {k ■ drf , (4.136) 

where the first term on the right-hand side is to be interpreted as in fiat space-time, with 
x^ = {x,y,z,t). Therefore u and k parametrize the deviation from flat space-time. 

Having introduced this new parametrization, the physical singularity — 0, equivalent 
to f = and 9 = tt/2, is now seen to be mapped to the ring 



85 



a;2 + y2 = z = 0, (4.137) 

with the solution 

X = — asiiK^, y = a cos if. (4.138) 

A consequence of this structure is, that in contrast to the Schwarzschild geometry, geodesic 
crossing the horizon do not necessarily encounter the singularity, but can pass it by. 



4.10 Black- holes and thermodynamics 

When a test particle falls into the ergosphere of a black hole, the constant of motion 
C,f P = £ becomes a momentum component, rather than an energy. This provides a means 
of extracting energy from a spinning black hole, as first proposed by Penrose [Q: if a 



particle enters the ergosphere with energy e, this becomes a momentum; during its stay in 
the ergosphere the momentum can be changed, for example by a decay of the particle into 
two or more fragments, in such a way that the final momentum of one of the fragments 
is larger than the original momentum. If this fragment would again escape, which is 
allowed as the ergosphere is located outside the horizon, it leaves with a different energy 
e' corresponding to its changed momentum. This energy is now larger than the energy e 
with which the original particle entered. By this process it is possible to extract energy 
from a black hole, even in classical physics. 

It would seem that the Penrose process violates the conservation of energy, if there would 
be no limit to the total amount of energy that can be extracted this way. However, it can 
be established that such a violation is not possible and that the process is self- limiting: 
emission of energy is always accompanied by a decrease in the angular momentum (and/or 
electric charge) of the black hole in such a way, that after extraction of a finite amount of 
energy the black hole becomes spherical. After that there is no longer an ergosphere and 
no energy extraction by the Penrose process is possible. 

We will now show how this comes about for the simplest case of an electrically neutral 
black hole of the Kerr type. In this case the location of the horizon is at 



th = rn + V ra^ — a? = m + — V — . (4.139) 

m 

Then the angular momentum of test particles moving tangentially to the horizon is 



VLh = , " , = — ] (4.140) 

This quantity can be interpreted in terms of the area of the spherical surface formed by 
horizon of the black hole: 



Ah = I I dOdip^/gM^ = 4vr (f% + a^) = (4.141) 

J Jr=ru ^ ' \Ih 
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Eq. ( [4.140| ) thus provides an expression for the area of the horizon in terms of the mass and 



the angular momentum of the black hole 



— = + - J2. (4.142) 

Now we observe that 

(i) the area Ah increases with increasing mass; 

(ii) at fixed area Afj the mass reaches a minimum for J = 0; 

(Hi) by absorbing infalling test particles (including Penrose- type processes) the black-hole 
area never decreases. 



The first two statements are obvious from eq. ( [4.1421) . The last one follows if we consider 
the effect of changes in the mass and angular momentum of the black hole. We assume 
that the absorption of a test particle with energy e and angular momentum i changes a 
Kerr black hole of total mass m and angular momentum J to another Kerr black hole of 
mass m' = m + e and angular momentum J' = J + i. This assumption is justified by the 
uniqueness of the Kerr solution for stationary rotating black holes. 

Now under arbitrary changes {6m, 6 J) the corresponding change in the area is 



, m 



5 Ah 2m5m / , , — —\ J 5. J 

(4.143) 



+ - J2) 
((5m — Qh^J) ■ 



^Vr - J2 V J _ J2 

J 



^H^m^ - J2 

For an infalling test particle only subject to the gravitational force of the black hole we 
have derived the inequality ( |4.109| ). Moreover, on the horizon such a particle has been 
shown to have the angular velocity VLh w.r.t. a stationary observer at infinity, independent 
of the value of Hence with 5m = e and 6J = i we obtain 

^ = , (e - nni) > 0. (4.144) 

Thus we can indeed conclude that by such adiabatic processes, adding infinitesimal amounts 
of energy and angular momentum at a time, the area of the black hole never decreases. 
Casting the relation ( |4.143| ) in the form 



Sm = —5 Ah + VLh6J, (4.145) 
Stt 

where k is given by: 

An I 

K = — — Vm^ - J2, (4.146) 
mAn 

it bears a remarkable resemblance to the laws of thermodynamics which states that changes 
of state described by equilibrium thermodynamics are such that for isolated systems vari- 
ations in energy U, entropy 5* and angular momentum J are related by 
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SU = T6S + Q6J, (4.147) 

where T is the temperature and fl the total angular momentum, and moreover that in 
such processes one always has 6S > 0. The important step needed in a full identification 
of these relations is to associate the area Ah with the entropy and the quantity k, called 
the surface gravity, with the temperature: 

5S = cSAh, T=^, (4.148) 

where c is an unknown constant of proportionality. Now if black holes have a finite tem- 
perature, they should emit thermal radiation. Using a semi-classical approximation, it was 
shown by Hawking ||4^, that quantum fiuctuations near the horizon of a black hole indeed 



lead to the emission of particles with a thermal spectrum. This is true irrespective of the 
existence of an ergosphere, therefore it covers the case of vanishing angular momentum as 
well. The temperature of the radiation Hawking found was 

Th = (4.149) 
This fixes the unknown constant above to the value 

c = ^. (4.150) 

It seems therefore, that the thermodynamical description of black holes is more than a 
mere analogy. The issue of interpreting the entropy 

S = ^ + So, (4.151) 

where 5*0 is an unknown additive constant, in statistical terms by counting microstates is 
hotly debated and has recently received a lot of impetus from developments in superstring 
theory. 

The extension of these results to the case of charged black holes is straightforward. The 
relation between horizon area Ah and angular velocity Qh remains the as in ( [4.141| ), but 



in terms of the macroscopic observables (m, J, e) associated with the black hole the area is 
now 

-Ji = + Vm^ - J2 - m2e2. (4.152) 

The changes in these observable quantities by absorbing (possibly charged) test particles 
are then related by 

Sm = —5 Ah + ^h5J + $H^e, (4.153) 
87r 

where the electric potential at the horizon is 
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= (4.164) 



and the constant k is now defined as 



47r 



K = — — Vm^ - J2 - m^e^. (4.155) 

Again, by sending test particles into the black hole the area can never decrease: 5 Ah > 0. 
In this way the thermodynamical treatment is generalized to include processes in which 
the electric charge changes. 

We conclude this discussion with establishing the gravitational interpretation of the con- 
stant K, associated with the temperature in the thermodynamical description. This quan- 
tity is known as the surface gravity, because in the spherically symmetric case it represents 
the (outward) radial acceleration necessary to keep a test particle hovering on the horizon 
without actually falling into the black hole. 

For the case of the Schwarzschild black hole this is easily established. Define the 
covariant acceleration as 

< = (4.156) 

As the four velocity of a test particle has negative unit norm: = —1, it follows that in 
the four-dimensional sense the acceleration is orthogonal to the velocity: 

ri-u ^ 0. (4.157) 

Also, a particle at rest in the frame of an asymptotic inertial observer at infinity, has a 
velocity four- vector of the form 

= (-^,0,0,o]. (4.158) 

Combining these two equations and taking into account the spherical symmetry, a particle 
at rest must have a purely radial acceleration given by 

- 1 Tn 

= -^dr-gtt = ^. (4.159) 

Remarkably, we recover Newtons law! It follows in particular that on the horizon 

r,^^± = K. (4.160) 
This quantity is actually a proper scalar invariant, as can be seen from writing it as 

« = V\Fa] , (4.161) 



rH 
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with T] the invariant magnitude of the four acceleration: 

V = ^Jgi^u^■ (4.162) 

For the case of the spinning black holes the proper generalization of this result gives the 
surface gravity as in eg. ( |4.155D . 
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